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An Investigation of the Physical Properties of Thin Films used 
in the Reduction of Surface Reflection * 


By. ac WEIS 
The Mullard Electronic Research Laboratory, Salfords, Horley 


MS. received 6th January 1948 


ABSTRACT. A simple mathematical treatment of reflection from thin films is extended’ 
to cover the properties of films likely to be met with in practice. In particular, reflection 
from multiple films of more than one substance is considered, and the reflection of polarized 
light from simple films. A photometer used for testing these results is described, and the 
different types of cryolite films found are enumerated. 

Evidence is given suggesting that polished or aged cryolite films behave as if they 
consisted of a lower stratum of cryolite surmounted by a secondary layer of higher refractive 
index. Some of the properties of this layer are described. 


§1. INTRODUCTION 


HE published theoretical treatment of evaporated cryolite films used for 
| reducing surface reflection from glass has not been extended far enough 
to explain observed results. A practical analysis based on the method 
used by Greenland (1943) has, therefore, been developed to explain the behaviour 
of films illuminated by polarized light, films with scattering centres within their 
bulk or irregularities on their surfaces, films of high dispersive power, and films 
consisting of layers of different refractive indices. Some of these latter conditions 
have now been treated by others (Mooney 1945, Vasitek 1947); the emphasis in 
the present paper is on a practical approach. 

A simple reflection photometer was used to compare the reflectivity of actual 
coated surfaces with that of normal untreated glass, over eight relatively narrow 
bands in the visible spectrum selected by colour filters. In many instances there 
were departures from the performance expected from simple theory which could 
only be explained by a more extended treatment; in particular, measurements 
suggested in some cases the production, either naturally or accidentally, of a 
secondary layer with relatively high refractive index lying on the original simple 
layer, and probably penetrating some way into it. ‘The presence of this Jayer was 
found to be of importance in producing an evaporated film which was of a stable and 
tough character. 


* Part of Thesis submitted to London University for Ph.D. Degree examination, 1947. 
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§2, MATHEMATICAL TREATMENT 
(i) General Method used 


In the mathematical treatment which was developed, only the condition of the 
reflected light was considered, as the construction of a transmission photometer of 
the necessary accuracy to check the behaviour of practical films was felt to be 
difficult. The analysis was based on the summation of multiple reflections, which 
shows that the intensity of light reflected depends simply on the effective reflection 
coefficients of the faces of a film having due regard to the path differences and 
transmission losses occurring within the body of the film. 

(ii) Single Films 

In the case of single films with no dispersion it can be shown that the intensity 
of monochromatic light J reflected from a given film varies with reciprocal wave- 
length or wave number 1/A, giving a curve of the form J~A+Bcosk/A where 
A and B are functions of the refractive indices present in the system, with A>B, 
and & a function of the optical thickness of the film. This simplified equation, 
due to the linear relationship between the wave number and the phase delay 
introduced when a wave traverses a film of constant thickness, was found to 
be theoretically accurate to about 5° in practical instances. 

By plotting experimental results in this way, recognizable fragments of cosine 
curves were usually produced from which the optical constants of the system 
could be deduced, although the measurements had only covered a relatively 
restricted waveband. 

(iii) Effect of Dispersion 

The effect of dispersion within a layer on a glass surface is two-fold: the 
difference in the optical thickness for different wavelengths gives a non-linear 
variation of phase with reciprocal wavelength and the reflecting power of the faces 
of the layer change with alterations of wavelength. 

A possible result of this is illustrated in Figure 1. Here at some wavelength A,, 
areflected amplitude given by X might be observed where X is a point on the curve 
(i). Atsome other wavelength A, the refractive index of the film might be higher so 
that a point such as Y lying on the new curve (ii) would be observed. Between A, 
and 4, intermediate values would be found, so that the curve (iii) is the effective 
reflectivity curve of the film. ‘The amplitude reflected could never actually 
exceed the amplitude reflected from a plain uncoated surface unless the film had 
such a high dispersion that the refractive index at some wavelength was greater 
than that for the glass base, but it is easy to see that incorrect maxima could be 
found if limited observations on curve (iii) were extrapolated without due care. 

If the dispersion of the glass is to be accounted for, further modifications will 
be produced, but only to the reflecting power of the second surface of the layer. 
‘This would raise or lower the minima of the respective curves, but it would only 
alter their reciprocal wavelength values by an amount of a second order. 


(iv) Multiple Films 


It was established that in the case of single films it was unlikely that a coated 
surface would be produced naturally which would reflect less light than an un- 
coated surface at some wavelengths and more light at others. Since such con- 
ditions were frequently found the possible existence of a secondary layer of a 
differing refractive index lying on or within a simple layer was investigated. 
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It was found as a general rule that if the refractive index steps of such a multiple 
layer system (even of more than three interfaces) were not all of the same sign 
{in going, say, from air to glass), then under some conditions it was possible that 
more light would be reflected than if only one step of the same overall magnitude 
existed. If all the steps were of the same sign, then the maximum possible 
reflection would not exceed that from a single step of the same overall magnitude. 


(v) Reflection of Polarized Light 


The reflection of polarized light to be expected from thin films was also studied. 
It was first assumed that no birefringent effect would make itself evident, in view 
of the thinness ot the films and the weakness of the birefringence of the materials 
used. As only parallel-sided films were considered, it was possible to consider any 
polarized incident light as made up of two quite independent components polarized 
in and perpendicular to the plane of incidence, for these two components would not 
at any time interfere with one another. ‘These two components were therefore 
investigated separately. 


(a) Light polarized in the plane of incidence 

The interference conditions for light falling at an angle on a single coated 
surface are similar to those for normal incidence except that the angles of incidence, 
and therefore the effective reflecting properties, differ at the two faces of the film. 


At the same time, the optical path within the film varies with the angle of incidence 
of the light. 


: 


RELATIVE REFL. AMPL. 


a (6 


Figure 1. Figure 2. 


For light polarized in the plane of incidence (magnetic vector) the amplitudes 
reflected at each face increase in a regular manner with increase of angle of 
incidence, but due to refraction within the film that from the upper surface will 
always increase at a faster rate than that from the lower face. ‘This results in a 
change in the relative intensities of the interfering components and thus in a change 
of saturation of the interference colour produced. 


(b) Light polarized perpendicular to the plane of incidence 
When the light is polarized so that the magnetic vector 1s perpendicular to the 
plane of incidence, the reflection obtained at each face of the film is a more complex 
function of the angle of incidence and includes the Brewster condition when zero 
reflection is obtained. Figure 2 shows the relative amplitudes of light reflected 
6-2 
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from the upper a, and the lower c, faces of a typical film with an index greater than 
the square root of the index of the glass on which it Jies, a condition almost always — 
met with in practice. These values of a, and c, are plotted as functions of «, the 
angle of incidence of the light on the upper face of the film. It is seen that a, 
decreases from its value at 0° to zero at Y, corresponding to the Brewster angle, 
then rises again but with an added phase difference of 7 radians. On the other — 
hand, due to refraction within the film, c, decreases very much more slowly and 
usually only just reaches the Brewster condition when « is near 90°. The optical 
path through the film will increase in the same manner as for light polarized in the 
plane of incidence. 

At X and Z, a, has the same value as c,, so that at these angles of incidence the 
correct interference conditions exist for the complete extinction of some wave- 
lengths in the reflected light. Thus at either point, the reflected light should have 
a saturated colour. 

The precise hue is determined at X by the path delay in traversing the film at the 
appropriate angle of refraction, while at Z the hue is determined by the path delay in 
traversing the film plus the additional phase retardation of 7 radians incurred on 
reflection at an angle of incidence greater than the Brewster angle. 

At Y, a, is zero, no interference is produced and no coloration of reflected 
light occurs. 


§3. PHOTOMETER 


The photometer used to measure the reflecting properties of thin films de- 
posited on glass was built almost entirely of Meccano parts; it is illustrated in 
Figure 3 (a) and its functional parts shown in Figure 3 (8). 

The image of the lamp L was thrown on to the surface to be tested 5, at nearly 
normal incidence (9°) by the converging lens C via a colour filter. A series of 
filters was used having fairly narrow pass bands distributed throughout the visible 
spectrum (Ilford ‘‘ Spectrum ”’ filters) so that substantially independent measure- 
ments could be carried out at differing wavelengths. 

A photocell P was arranged to lie within the limits of the reflected cone of light 
at such a position that quite wide tolerances in the position of the specimen were 
posssible ; this simplified the changing of the specimens. For rapid test a rotating 
table was used on which three coated specimens and a plain comparator could be 
mounted. 

Since the light to be measured was extremely faint the photocell used was a 
Cinema-Television multiplying cell, which was arranged to multiply the photo- 
current emitted by the cathode by rather more than 104 by the process of secondary 
emission. A sensitive mirror galvanometer (5 ya. full-scale deflection) in the final 
collecting electrode circuit indicated the intensity of the light. The Cinema- 
Television caesiated-antimony cell surface, insensitive to infra-red, was used, since 
colour filters were found to be transparent to the near infra-red. 

From the published data on the eight colour filters, the photocell and the light 
source used, it was possible to associate with each filter a wavelength of maximum 
transmission, which was taken as characteristic for the measurement. In spite 
of the difficulties introduced by the finite pass bands of the filters, much greater 
sensitivity was thus obtained than would have been possible with a monochromator. 

The reflection from a specimen was compared with the reflection obtained 
under identical conditions from an uncoated surface of a similar glass block. The 
unwanted reflection from the second glass—air surface of each block was removed by 
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Figure 3 (a). Photometer 


Figure 3 (0). 
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were variable factors, differing from specimen to specimen, but although in some 
cases little or no displacement was produced, the rise in the minimum was always. 
found to exist, and often lay between 2% and 5%. 


(ii) Reflection of Polarized Light 

With a simple polarizing spectrometer it was possible to confirm the theoretical 
reflecting properties of thin films at angles of incidence other than normal. The 
usual prism was replaced by a specimen arranged to reflect light from the collimator 
down the telescope tube; the appropriate angles could then be read off in the usual 
way. No trace of elliptical polarization was ever detected, the observable limit 
being about 1%, so that it was considered that the elementary treatment was 
justified. 

(a) Light polarized in the plane of incidence. When films were illuminated 
with light polarized in the plane of incidence (magnetic vector) the hue of the light 
reflected changed slowly with increasing angle of incidence as was expected from 
the increasing optical path within the film. The colour saturation, however, 
with all types of films decreased slightly with increasing angles of incidence, 
indicating that the index of refraction of the film was always rather higher than the 
square root of the glass base. 

(b) Light polarized perpendicular to the plane of incidence. ‘The effects observed. 
when films were illuminated with light polarized perpendicular to the plane of 
incidence (magnetic vector) were also as expected. Normal films were found to. 
show a hue which changed slightly as the angle of incidence was first increased,. 
and which had an almost velvet-like appearance when the critical state was 
approached. Colour saturation occurred usually at some angle of incidence near 
to 40°, after which the saturation faded rapidly and the film reflected nearly white 
light over a large range of angles. Nearer to grazing incidence it was sometimes. 
possible to see the approach of the second saturation condition, but this was not: 
always obvious and was never observed at any definite point. 

When some colour could be seen, however, the hue of the light reflected when _ 
the plane of polarization was in the plane of incidence was complementary to that 
observed when the plane of polarization was perpendicular to this, at the same 
angle of incidence. This was caused by the extra phase change which occurred in 
the latter case when the angle of incidence of the light on the air surface of the film 
was greater than the Brewster angle. 

It was found that the angle of incidence « at which colour saturation occurred. 
could be measured to an accuracy of about 1°, though with some films, colour 
filters in the eyepiece of the spectrometer were necessary to attain this figure. 
From these measurements the refractive index of a film could be calculated assum= 
ing only the index of the glass base. A comparison of results obtained by this. 
method and by the photometer method described in (i) is given in two cases : 
for a soft film the photometer method gave minimum reflectivity 1°% and refractive 
index 1-285 while the polarization method gave a refractive index 1-275 from a 
value of «=30°; for a lightly polished film the two values were 1-367 and 1-355, 
from miminum reflectivity 14% and « =46°. 

In spite of the good correspondence, it was found from later work that it was not 
always possible to assume a simple layer of uniform refractive index in practice, 
so that the convenience of this single observation method was lost. In general, 
however, a high value of « was found to indicate a high value of refractive index. 
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(ui) Effect of Polish 

(a) General. It was soon found in experiments on thin films that the effects 
of polishing were profound. ‘The most striking observation made was that as a 
fresh, unpolished surface was rubbed, a characteristic squeak was produced until 
all the surface had been disturbed, when no amount of further rubbing, even after 
a period of months, would produce the squeak again. Tests indicated that the 
mechanical resistance to the polishing material was much higher before the squeak 
than after it and also that dusting or rubbing lightly with a paint-brush or cotton- 
wool did not cause the effect, or prevent it being found afterwards. It was found 
in addition that a polished surface was much more robust than an unpolished one. 
Thus, cotton-wool rubbed briskly over a fresh surface would partially remove it, 
while producing the characteristic squeak, but the remaining film, being polished, 
would then be very stable and resist further rubbing without even scratching. 

These effects were confirmed with many polishing materials, but while silk 
and cotton were both found to polish without substantial initial scratching, 
chamois-leather was even better in this respect and was therefore more often used 
when photometric measurements were made. 

To check these findings, a simple scleroscope was made. ‘This instrument was. 
not used to measure the abrasion resistance of the films, as did that due to Townsley 
(1945), but to produce, under carefully controlled conditions, a single scratch 
which could be examined microscopically. The results suggested that while on 
unpolished films fine continuous scatches were produced, when the same films 
were polished the scratching point tended to skid over the surface making a 
regularly spaced series of indentations into the film, doubtless due to point chatter. 

These mechanical properties could have been explained by assuming that the 
squeak corresponded to the cleavage and pulverization of the oriented crystallites. 
found by electron diffraction in an unpolished layer, and always absent or masked 
in a polished specimen (Dell, in preparation) but the optical properties produced at 
the same time showed that this explanation was not sufficiently detailed. 

When the reflectivity results given by the photometer were extrapolated by 
fitting cosine curves to them, it was found that unpolished films usually gave curves. 
which had maxima corresponding to a reflectivity of rather more than 100%. A 
typical result is illustrated in Figure 5, curve A, where a maximum of 116% was 
obtained. This value is not altogether surprising when the uncertainty of such an 
extended extrapolation by an approximate method is considered. 

Whena similar attempt was made to fit a curve to the reflectivity results ahiinerl 
from the same specimen after polishing (curve B) a maximum of 126°% was found. 
This was outside the experimental error expected, and a confirmatory check was 
carried out on a film (Figure 6) of such thickness that both a maximum (second 
order) and a minimum appeared within the visible spectrum. ‘This particular 
specimen was deposited on glass of refractive index 1-61, and from the minimum of 
the curve found when the film was unpolished (curve A) it can be seen that the film 
index was about 1-34, while the maximum reflection corresponded to about 110%. 
As soon as the film was polished, however, the maximum rose to 150% (curve B) 
while the minimum shifted in the usual manner, indicating that the effect suspected 
in thin films certainly had a counterpart in a thick film. 

In the experiments on polishing it was found that there was a great danger of 
contaminating the Pe suniaee 2y grease, when the polishing action tended to 
produce a marked bluish “smear” on a normal specimen which could only be 
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partially removed by brisk polishing afterwards with cotton-wool. ‘Thin rubber 
sheeting wrapped round a cotton knob was found to be satisfactory except that a 
smear of greyish matter was deposited on the film, but this was removable by 
subsequent rubbing with cotton-wool. This left the film quite hard and polished 
to a high degree with no visible scratching. ‘Typical reflectivity curves were 
found to be very similar to those produced by polishing with chamois-leather, 
except that the rise in the minimum reflectivity was usually rather greater, so that 
16% to 17% of that of plain glass was common. ‘The extrapolated maximum, 
however, was rarely greater than about 110%. 

(b) Existence of secondary layer. It has been shown that if simple dielectrics 
alone are considered, reflectivity curves with maxima greater than 100° are only 
possible when some double layer is present in which the steps of refractive index 
are not all of the same sign. Thus it was likely that polishing a film produced 
either a zone of abnormally high refractive index at the surface or a zone of ab- 
normally low refractive index near the glass, or possibly both these effects. 
Experiments were therefore carried out to determine whether such changes were 
present. 

(c) Special tests on polished specimens. 'The experiments on hardness suggested 
that of the two possibilities the disturbance of the upper surface was the more 
probable, and from the effect of grease on the polishing material it was thought 
possible that grease solvents or possibly some etching technique might be found 
which would remove the suspected secondary layer either partially or completely ; 
little success was obtained along these lines, but it was found that when a polished 
film was flooded with a solution of collodion in acetone, it was possible to strip the 
collodion layer off the film when it had dried out, and when this was done, some at 
least of the properties of the unpolished specimen were restored. 

The scleroscope indicated that stripped areas of polished films were scratched 
exactly as new, unpolished films; zones which had been stripped stood out clearly 
as they reflected noticeably different and more saturated colours than the un- 
stripped polished zones. Only in the fact that fresh polishing failed to produce 
any squeak did they differ from fresh specimens. 

Measurements illustrating the history of a typical specimen repeatedly 
polished and stripped in this way are given in the table which shows the relative 
reflected intensity J, and the reciprocal wavelength of the experimentally deter- 
mined minima and the extrapolated maxima (the value of J, for plain glass is 
100). Curves A~G were drawn at various stages and the results obtained from 
them are summarized in Figure 7. 

The newly deposited film gave curve A and immediately after this had been 
polished by the rubber process, the curve B was produced, and 36 hours later 
curve C. Further tests indicated that no further changes were occurring. 

Three years after this the specimen was cleaned vigorously with a cotton 
handkerchief and curve D was plotted, and from this it seems clear that a sub- 
sidiary layer was then present. ‘This film was then stripped (curve E), re-polished 
with cotton (curve F), and stripped again (curve G). 

This complete cycle of changes is summarized in Figure 7 which shows the 
reciprocal wavelength of the minimum and the actual value of that minimum 
after each significant action. As can be seen, each period of polishing raised the 
minimum and pushed it to a smaller reciprocal wavelength, implying an increase 
of the total optical thickness, while each stripping lowered the minimum and 
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‘increased the reciprocal wavelength at which it occurred, ane a reduction of 
-optical thickness. 

_ These results were not exceptional though in a few cases it was not found easy 
to strip a film evenly, possibly due to faulty deposition or poor attachment to the 
glass, but, judged visually, the results were the same even in these cases. 

Another unexpected phenomenon almost certainly of an allied nature was 
found in the case of fairly thick specimens which had originally been of the low 
index type. One of these had a maximum of almost exactly 100% when first made 
and a minimum of about 5%, showing that the index was then 1-32, but after 48 
hours it was found that a curve of the type A in Figure 8 was produced; asis seen 
this had a maximum of more than 200%, the highest Figure ever found. In order 
‘to determine the crystalline state of this untouched surface, the specimen was 
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Figure 7. 
Curve Minimum (experimental) Maximum (extrapolated) 

Iz po Tp po 
A 10-6 2:07 sales} 4-14 
B 17 1-84 114 3-68 
Cc Gis 1:85 — — 
D 18 19/7 128 3°54 
E 11 1-85 109 3°7 
F 13e5 1-85 119-5 Bez 
G 12°5 1-88 112°5 3-76 


‘placed in an electron diffraction camera, but under the influence of the vacuum 
-conditions or of the low velocity (150 v.) subsidiary electron stream, the surface 
“peeled” quite suddenly, forming a fine silvery dust of curled flakes of film lying 
-on the glass. When these were dusted off with a brush, it was found that a portion 
-of the film remained behind, leaving a perfectly even pale brownish-coloured 
layer, which gave the reflectivity curve B in Figure 8. 
Little can be deduced from this latter curve, but it is possible that it represents 
‘the lower portion of the multiple film whose existence was implied by curve A 
_and which might have remained after the sudden removal of the upper component. 
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The only other specimen which stripped itself in this way did not leave a surface 
suitable for photometric work, but judged by eye it was rather thicker than the one 
illustrated. 


(d) Inferences. 


Probable nature of layer. All these experiments indicated that it was quite 
likely that the changes found in evaporated films when they were polished or 
allowed to age could be explained if it were assumed that both these processes 
tended to produce on the free surface of the films secondary layers of a relatively 
high refractive index. From the restricted range of wavelengths observed by the 
photometer it was difficult to gain much information about the optical properties. 
of the secondary layer and it was even uncertain whether it lay on or within the 
boundaries of the original film, but some suggestions may be made as to its. 
probable properties. 


150 


100 


50 


Relative reflected intensity (plain glass—100) 


1:5 1:75 2:0 2°25 
Reciprocal wavelength (u-1) 


Figure 8. 


In Figure 6 (B) a maximum reflectivity of 150% is shown; the minimum. 
refractive index of the subsidiary layer necessary to produce this reflectivity can 
be found from the theoretical analysis. If it is assumed that the refractive index of 
the cryolite was 1-34, even after polishing, and that the phase relationships were 
such as to produce reinforcement by all components at the wavelength of maximum 
reflection, then it can be shown that the refractive index of the extra layer must be 
1-4, If, as is probable, the phase relationships were not so favourable, a higher 
index would be necessary. Similarly, in Figure 8 (A) the maximum of about 
200% indicates an index of at least 1-43 for the extra layer, if such a single extra. 
layer exists. 

Assuming that such layers were possible, it was necessary to find whether the 
curves as a whole were of the form that might be produced by such means. 
Theoretical curves were therefore drawn showing the effect on the reflectivity 
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curve when different thicknesses of material of a relatively high refractive index. 
lay on the surface of a normal film. In Figure 9 are shown approximate curves 
which give the general trend of the results, and agree with the exact formula to 
an accuracy of better than 5%. 


In Figure 9 it is assumed that a film of cryolite of 0-1 thickness and refractive. 


Relative reflected intensity (plain glass—100) 


0 100° 200° 300° 400° 
Phase delay in lower layer (83). 
iy— le (air), nN, =1-36 (cryolite)—O-1p thick, 
My = 1-5—various thicknesses. N,=1°58 (glass). 
Curve i 2 83 4 5 6 
Thickness of m.(a.) 0 50 150 250 500 907 
Figure 9, 


index 1-36 lying on glass of refractive index 1-58 is covered by a layer of material 
of index 1-5. Reflectivity curves are plotted as a function of phase for thicknesses 
of the subsidiary layer between 0 and about 0-1. The portions of these curves 
which would be observed by the photometer lie between the points X and Y on 
the phase scale. For other thicknesses sections of the curves observed would lie 
at other points along the phase axis, scaling the actual thickness of both layers in the 
same manner. 

It can be seen from these curves that the rise in the maximum and minimum 
and the change in their positions as the subsidiary layer increases in thickness are all 
of the order of magnitude that were found in practice; in particular, the extreme 
sensitivity of the reflectivity curve is shown when even as little as 50 a. of any extra 
material is added to a surface. 

An experimental fact not yet satisfactorily explained is that often when polish- 
ing was carried out (Figure 5) the shift in the minimum found was less and the rise 
in the minimum more than would be indicated by a simple process such as has been 
visualized here. While this explanation might be partly right, therefore, it does. 
not entirely cover the facts. On the other hand, the marked departure of the 
curves shown in Figure 9 from a simple cosine law explains the rather low values. 
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that were sometimes found for the maximum reflectivity of polished specimens 
when these were calculated by extrapolation assuming such a simple law. 


Possible causes of layer. It is not difficult to imagine how such additional layers. 


could come to lie on the surface of the cryolite, for many animal and vegetable fats 


and oils have refractive indices between 1-47 and 1-48, and the manner in which 


such oils can creep over surfaces is wellknown. Thus even the vaseline (1 = 1-49) 
used to reduce unwanted reflection, or the pump oil from the diffusion pumps 
(n= 1-48) could be the origin of the ageing effects, while the effects of polishing 
might well depend on rubbing some grease such as cotton seed oil (n = 1-475) into 
the upper layers of a film. Chamois-leather in particular (Woodroffe 1941) 
consists of animal skins tanned with unsaturated fish oils (n= 1-48) and these, 
though often partly oxidized by exposure to the air, could almost certainly form a 
very close link with the cryolite molecules in the surface of the film, giving a 
particularly stable layer. 

These facts would explain, too, the very great difference in mechanical pro- 
perties between the polished and unpolished specimens, and no doubt a smooth, 
grease-protected surface would be in a much more permanent and stable state than 
newly evaporated cryolite. This improved stability of polished films was very 
marked in the case of moisture-provoked changes, as it was found that well- 
polished specimens were often completely unaffected by water vapour from any 
source. 

It is not at once obvious in just what way rubber polishing operates but there is 
no doubt that the deposition of any form of contamination of a relatively high 
refractive index would inevitably change the reflection conditions in the manner 
shown. It is probable, therefore, that even after cleaning off the grey smudge 
produced by such polishing, a thin protective layer of rubber remained behind. 
Too few experiments have been made on rubber-polished films to obtain much 
information about them, but doubtless the effect was similar in nature if not in 
detail to polishing with other materials. 

It was found necessary to assume that in many cases the additional layer had an 
appreciable thickness and it is strange that even layers of the order of 250 a. or 
more were so well-attached to the cryolite that repeated heavy rubbing with other 
polishing agents did not always change them. The self-stripped specimens 
suggest, however, that the high-index layer may sometimes penetrate some 
distance into the original cryolite film, since this would explain the rupture of the 
specimens under the action of the discharger, a known degreasing process. It is 
just possible that oil or grease could soak into the upper regions of a film, probably 
increasing the mechanical dimensions as it did so, but making a relatively strong 
mass rather than a weak layer of mobile oil on a hard surface. 

In addition to this, if such soaking existed it is probable that the simple 
condition suggested of two films each of uniform optical density would rarely 
occur in practice. It would be far more likely that the different materials would 
merge into each other so that more complicated conditions of reflection would 
exist. These would probably explain the difficulties that have been found, but 
their theoretical treatment is not simple. 
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ABSTRACT. ‘The single reflection echelon grating has the disadvantage that it is 
impossible to interpret the fringe pattern without further information about the correct 
order allocation. ‘The present methods of obtaining this information are briefly reviewed ; 
they are progressively less useful as the wavelength decreases into the ultra-violet. 'The 
properties of the grating are investigated theoretically when the light is incident at relatively 
large angles in two special cases, and the performance of the instrument in these positions 
compared with that of its normal use. It is shown that the required order allocation could 
be obtained using a single reflection echelon grating by measuring the fractional order 
separation with the grating in two positions. This would considerably increase the value - 
of the instrument. 


§1. INTRODUCTION 


HE importance of the reflection echelon grating in the field of high resolution 
spectroscopy is due partly to the fact that it has a good uniform performance 
over a wide range of wavelength. But, more especially, it can be used for 

the determination of absolute wavelengths and of small wavelength differences in 
the ultra-violet where all other high resolution instruments cannot be used: for 
various reasons. ‘The only exception to this is the vacuum concave grating whose 
resolving power is inferior to that of the other instruments. 

However, when the reflection echelon grating is used in the ordinary manner 
with crossed auxiliary dispersion, it is impossible to interpret unambiguously the 
fringe pattern, i.e. to give the fringes their correct order numbers, the order of a 
fringe being the integral number of wavelengths in the optical path difference of 
two consecutive interfering rays. Also it is impossible to detect the coincidence 
of the fringe systems of different components when this occurs. ‘Thus until 
additional information is obtained, the true character of the spectral complex 
under consideration cannot be deduced. 

It is proposed here to extend the theory of the reflection echelon grating and 
to show how, using one such instrument, it would be possible to make suitable 
observations which allow the correct interpretation of the fringe system to be made.. 
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§2. REVIEW OF THE PRESENT METHODS 
OF FRINGE INTERPRETATION 


We show in Figure 1 the type of intensity distribution we should get with a 
simple spectral complex consisting of two components of unequal intensity. ‘The 
intensity envelopes are shown for reference. The most general order and wave- 
length allocations are given below the diagram. 


Wavelength A Ardr a A+dr 
Fringe Order =m m’ (m+!) (n't |) 


Figure 1. Intensity distribution for a two component complex. 


The fractional order separation of the fringe of order m’ of the component 
(A +dA) from the fringe of order m of the component A is given by f=x/y, where x 
is the mean separation of two adjacent fringes of different components, and y that 
-of two fringes of the same component. 

If we use the fact that the product of the order and wavelength is a constant 
for any one angular position, we get the following relation 


dA = (tt —micef Ali 9 ee (1) 


The one plate alone gives no information about the value of (m—m’), and dX is 
-one of several possible values; (m-—~m’) will be referred to as the overlapping. 

The information needed to resolve this uncertainty may be obtained by two 
essentially different methods. 

Probably the more convenient method is to use another high resolution 
interferometer whose fringes are of considerably different order ofinterference. If 
the fractional order separation is measured, the correct order of the overlapping is 
easily deduced. ‘lhe instrument may be used separately or combined with the 
reflection echelon grating and its performance need not be as good. However, 
the limitations of the various interferometers restrict the use of this method to 
wavelengths greater than 2000a. We shall not consider the use of a second 
reflection echelon grating, as few laboratories have two such instruments available 
and we wish to investigate the methods of analysis which use only one of these 
instruments. 

Alternatively the reflection echelon grating and the auxiliary dispersive 
-apparatus may be mounted so that their dispersions are parallel. They may then 
be set so that their dispersions are assisting or in opposition and there is a small 
. change in the fractional order separation in the two positions. This shift can be 
used to deduce the value of the overlapping. The shift in the ultra-violet using 
a quartz spectrometer is easily detected; Williams (1948) observed a shift of some 
4% in the near ultra-violet using an E1 spectrometer. But the lower wavelength 
limit of use of a quartz spectrometer is in the region of 2000 .; for shorter wave- 
Jengths a diffraction grating must be used. There is no recorded instance of its 
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) use up to the present time, but a simple calculation gives the following formula for 
the change in fractional order separation in the two positions using a concave 
| grating in the Eagle mounting : 


Af ={Mns/[1 —(Mna/2)?]?}(m — m' +f )a/t, - 


/ where M™ is the order of the concave grating in which observations are made, 
} n the number of lines per centimetre of the grating, s the offset of the echelon 
| grating, and ¢ the echelon plate thickness. The substitution of typical values 
indicates that the method would be unreliable for wavelengths shorter than 1750 a. 


The parallel type of mounting, moreover, has the general defect that it can be 


used only when the number of spectral lines present is small and they are well 
» spaced. ‘Thus the method is applicable in a strictly limited number of cases. 


Hence it would appear that while the present methods of employing a single 
reflection echelon grating are satisfactory in the visible and near ultra-violet 


- spectral regions, they become progressively less useful as the wavelength decreases. 


We are therefore led to investigate the possibility of using a single reflection 


echelon grating in a suitable manner which would enable the correct order allo- 


cation to be made, and the true nature of the spectral complex derived. 


3 LHEORY OF LHE RERLECTION ECHELON GRATING 
The reflection echelon grating is usually mounted so that the light is incident 


| almost perpendicularly upon the plate offset, the incident and diffracted beams 
| making small angles with the normal. We now wish to investigate its properties 


when the light is incident from two other directions. For convenience of descrip- 
tion we shall refer the grating to rectilinear axes, the yz plane being parallel to the 
plate offsets with the z direction parallel to their long edges; the x direction is thus 


_ that of the offset normal. ‘Then, in the usual mounting, the light is incident in the 
_ .xy plane and makes a small angle with the x direction. 


(i) Light Incident almost in the xz Plane and at Large Angles of Incidence 


The general directions of the incident and diffracted beams are shown in 
Figure 2. The angles« and £ are the angles made by the incident and diffracted 


Figure 2. 


(a) Section in the xz plane. (b) Section in the xy plane. (Rays not in the 
plane.) 


rays respectively with the xz plane and are positive when on the same side of the 
normal. From the geometry of the figure the path difference between correspond- 
ing rays at adjacent offsets is 2t cos 6(cos « + cos 8) —s(sina + sin f). 

The general expression for the total intensity distribution from N regularly 
spaced offsets is [=J,(sin Nd/sin¢)?, where 2¢ is the phase difference between 
corresponding rays and J, represents the intensity distribution from one offset. 
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The intensity distribution of the light diffracted from one offset is given by the 
expression J,=J)(sin#/s)? where ~=n(sina+sinf)s/A and I) is a constant.. 
Thus in this case, we have J = J,(sin %/s)?(sin Nd/sin d)? where 


2h =27{2t(cos « +.cos B) cos 6 —s(sin« +sin f)}/A. 


In practice, the angles « and f will be very small so that the usual trigono-. 
metrical approximations (sina«=«, etc.) can be applied safely. Hence for the 
principal maxima 2¢ cos @—s(«+)=maA or defining an angle A by A=a+f, we 
get 2tcos@—sA=mdA; A represents the angular displacement from the rays 
obeying the laws of reflection. 

The following formulae for the grating are easily obtained by differentiating 
this expression and using the approximate relation 2fcos#=mA where appro- 
priate : 

Dispersion : variation of A with A at constant order number 0A/0A = 2t cos 8/sA. 

Dispersion of orders : variation of A with order of interference for one wave- 
length dA/dm=A/s. 

Wavelength range (A’A) between successive orders: (A’A)=A?/2tcos@. This 
is the maximum wavelength difference possible between two components so that 
no overlapping occurs. 

Resolving power: A/SA = Nm =2Nt cos 6/A, where 6) is the smallest wavelength 
separation measurable at a wavelength A. 

These formulae are very similar to those of normal incidence which are 
Dispersion OA/OA = 2t/sA Wavelength separation (A’A) =A?/2¢ 
Dispersion of orders 0A/dm=A/s Resolving power A/SA=2Nt/A. 

The graph of Figure 3 shows the variation of the performance of the reflection 
-echelon grating with the angle @. In each case the corresponding value at normal 
incidence has been taken as being unity. The broken line gives the variation of a 
quantity which we shall call the fractional change in path difference. This is 
given by (L— L’)/L, where L is the path difference at normal incidence and L’ that 
under the given conditions. ; 

It will be seen from the graph that there is a relatively small deterioration in. 
the performance of the instrument. One point of particular interest is that the 
dispersion 0A/0d is still linear, i.e. independent of A. 


(ii) Light Incident in the xy Plane at Large Angles of Incidence 


This is really the general case of the usual method of mounting as the light is 
incident in the same plane, but no restrictions are being made on the magnitude of 
the angle of incidence except that it shall be so large that all rays incident on the ~ 
offset, if reflected, would be incident upon the plate side (Figure 4). Expressed 
symbolically, this is tan~!@>s/t where s and ¢ are the plate offset and thickness. 
respectively. 

Considering any two corresponding rays at adjacent offsets, it is easily shown. 
that the path difference between them is given by ¢(cos@ + cose) —s(sin @ +sine) 
and is independent of the intermediate angle of diffraction y. 

The resultant intensity of illumination in the focal plane of the collecting lens in 
the direction ¢, is given as before by 


I=I,(sin N¢/sin d)?, 
where 2h = 27 {t(cos 0 + cos) — s(sin 8 + sine)}/A. 
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| The intensity distribution J, from the single element is not so easily found 
i exactly as in the case of the long reflecting strip. However, an elementary’ 
| consideration of the physical principles involved indicates that the distribution 
| will be similar to that of a reflecting strip of width 2s, the light being incident at an 
| angle 6 and diffracted at an angle e. The actual shape of the intensity envelope is 
| of little importance provided there are no unusual singularities present and that the 
| angular separation of the central minima is sufficient to contain two orders of each. 
| component. 

Thus the intensity will be given by J=J)(sin:/)?(sin Nd/sind)? where 
#=2n(sin@—sine)s/A. For the principal maxima of the diffraction pattern, we 
: have t(cos 6 + cose) —s(sin@ + sine) =ma. 


10 20 
Angle of Incidence (degrees) 


Figure 3. Variation of the relative performance 


with the angle of incidence. Figure 4. Section of the echelon grating in 


. Wavelength separation of orders 4’A. the «xy plane. 


. Dispersion of orders 04/0m. 
. Resolving power and dispersion 04/0A. 
. Fractional change in path difference. 


PRONFR 


If we define an angle A through the relation A =¢ —6, a consideration of the 
magnitudes involved shows that we may write sin « = sin (6 + A) =sin@ + Acos @ etc. 
where A represents the angular displacement from rays obeying the normal laws of 
reflection. ‘Then for the principal maxima we have 


2(t cos 6 —s sin @) — A(t sin @ +s cos #) = mA. 


The following formulae for the reflection echelon grating are easily obtained by 
differentiating this equation.and using the approximate relation 


mA = 2(t cos 8 —ssin 6) 
where appropriate : 
Dispersion, i.e. at constant order number: 


dA/O0A = 2(t cos 6 —s sin 6)/A(t cos 8 +s sin @). 


Dispersion of orders : 0A/Om=A/(tsin 6 +s cos 6). 

Wavelength range between successive orders : A’A =A?/2(t cos 6 —s sin 6). 

Resolving power : 4/8A = Nm =2N(t cos @ —s sin @)/A. 

The performance of the instrument is shown in Figure 5. As in the previous 
case the broken line gives the variation of the fractional changein path. The value 
10 has been assumed for the ratio of plate thickness to plate offset. However, it is 
only a relatively unimportant factor and the curve will be very similar for any other 
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value of the ratio. The graph indicates clearly that while there is no great deterior-_ . 
ation in the resolving power, the dispersion falls away very rapidly. his is likely a 
to be a very serious defect in view of the small dispersion of the instrument even in — 
its usual mounting. i 


‘ 
* 
j 
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§4. USE OF A SINGLE REFLECTION ECHELON GRATING 
FOR ORDER ALLOCATION Ms 

Without considering at this point the relative merits of the two possible 
methods of use investigated, we shall see how it would be possible to use either of 
them to decide the question of the correction order allocation. 


I-k , ie 
3 
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E 
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S 
: : “# 
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es : 
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Figure 5. Variation of the relative performance k -2 
with angle of incidence. = 
1. Wavelength separation of orders (4’)). % 
2. Resolving power. : f 
3. Dispersion of orders 04/dm. Figure 6. Graph showing the relation between 
4. Dispersion 04/0A. the measured fractions in the two positions. 


The reflection echelon grating is first used in the normal mounting and — 
accurate measurements are made on the various fringe separations. ‘Then when~ 
the question of the order allocation is decided, the wavelength distribution and 
separations of the components can be derived. 

The essence of the method proposed is to alter the path difference between rays — 
from adjacent offsets by having the light incident from one of the two directions 
previously investigated, and then from the new fringe separations we can deduce 
the correct order allocation. 

Considering a fringe distribution similar to that of Figure 1, and giving the same 
general order allocation, the measured fraction is related to the component wave- 
length separation by the expression 


dA = (m= PN Uti | tea. ie eee (2) 
as previously given. Now m, the order of interference, is connected with the 
path difference between corresponding rays by the relation mA=L, and thus 
equation (2) becomes 

dh=a(m—n tpl e 7 Sa eee (3) 

Now supposing that we use the echelon grating in such a manner that the path 
difference is altered to some smaller value L’. This will cause a change in the 
fringe pattern and the new measured fraction f’ will be given by 


BN = (10 JL ~ at rarer (4) 


where the n’s are the new orders of interference of the two components. 
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We shall consider the measured fraction. between two components in the 


| instance the component which we designated (A+ dA), and the nearest fringe of the 
) other component on the shorter path difference side. ‘This is by no means the only 
way in which we can define the fractional order separation and, in fact, the actual 
) definition adopted is of little consequence provided the same method is always used. 

_- Then from equations (3) and (4), we get 


fi=(m—m)L'/L-(n-n)4 fl [k,n (5) 


The relation between the two measured fractions given by equation (5) is best 
)shown graphically by plotting /’ against f for various values of (m—m’'). The 
» corresponding values of (m —n’) are easily deduced from a consideration of equation 
( (5) as we decrease the value of L’/L from unity. In the graph, k has been used to 
represent (L—L’)/L for convenience. ‘The graph is easily constructed for any 
particular value of L’/L by marking off on the axes the points indicated and 
drawing parallel lines through them with a gradient equal to L’/L. The dis- 
continuities which occur in the graph result from our definition of the fraction to 
_be measured and correspond to a change of unity in the value of (n—n’). At the 
discontinuity, the fringe systems of the two components coincide. Thus to get 
the correct order of the overlapping (m—m’), we have only to find from the graph 
the value which gives the correct value for the measured fraction f’ corresponding 
to the original fraction f. 

If the components are nearly of the same intensity, it is quite likely that the 
fraction measured in the second position will not be the fraction satisfying our 
definition but that given by f’=1-—/’. This, in practice, will lead to little trouble, 
if any, as we shall find that there is no fraction of this value on the graph correspond- 

ing to the original fraction and thus we should naturally try the fraction given by 
af"). 

_ The fractional aera of two different values of the overlapping in the 
second position is equal to (L—L’)/L so that if we can measure the fraction to 
an accuracy of about one quarter of this amount, we can always derive the correct 
order of overlapping. Thus if we can measure a fraction to one hundredth of an 
order, the fractional change in path difference can be made equal to a twentieth or 
so. A value of a tenth would probably be ory convenient as a fixed position in. 
which the second fraction is to be measured. 

In the rarer case when we have the coincidence of the fate patterns , of two, 
different components in the original position, this will be at once obvious in the 
second position by the appearance of another pair of fringes. ‘The large fractional 
separations of different values of the overlapping ensures that the components will 


now be clearly resolved. 


§5. DISCUSSION OF THE PROPOSED METHODS OF USE 
Any method of using the reflection echelon grating so that the path difference of 
corresponding rays is different in the two positions will be theoretically suitable for 
the allocation of order numbers. The position in which the fringe measurements 
must be made with high accuracy will be that in which the grating is used in the 
normal manner and the graphs of Figures 3 and 5 show clearly that this is the 
position of optimum performance. In the second position, no great accuracy is 
if the fractional change in path difference is equal to a tenth, the 
7-2 


required : 
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fractional fringe separation need only be measured to one-fortieth of an order. 
Therefore, a relatively large deterioration in the performance of the instrument is 
permissible. 


The first method of having reflections from the offset and from the plate edge — : 


has the advantage that the diffracted beam returns along the path of the incident — 
light so that the instrumental set-up can be very similar to that normally used, the 


only modification needed is some arrangement for tilting the echelon grating — | 


through the required angle. But the disadvantages are so great that this method 
will probably not be practicable. It is highly improbable that the plate edges are 
polished to the high standard required and adjusted perpendicularly to the plate 
offset and the cost of so doing would be very great. At the same time, the dis- 
persion falls away very rapidly with the increase in the value of the fractional — 
change in path difference, and in view of the small angular dispersion of orders of _ 
the instrument even in its most favourable position this becomes a prohibitive _ 
factor to the use of this mounting, especially as it is most likely to be needed in the ~ 
ultra-violet where the dispersion is least. 

The second method appears to be much more promising. ‘The performance 
of the instrument remains very good and the small decrease in dispersion is quite — 


tolerable. Further, the performance depends only on the surface finish and 


accuracy of the plate offsets and this must be good for the normal use. ‘The great 
disadvantage is that the diffracted light does not return along the direction of the — 
incident light and thus the ordinary type of mounting cannot be used. There 
appears to be no practicable arrangement of mirrors which will reflect the diffracted | 
beam along the direction of the incident beam. The grating would probably have 
to be mounted in a manner similar to that of the plane diffraction grating. This is 
unfortunate in that it involves another lens and a complete rearrangement of the 
auxiliary dispersive apparatus, yet as we can tolerate some deterioration in per- 
formance, the quality and adjustment of this apparatus need not be as good as 
that of the principal apparatus. 


§6. CONCLUSION 
The theoretical investigation of the properties of the reflection echelon grating 
shows that the performance remains relatively good as the angle of incidence in the 
plane parallel to the plate edges is increased. Measurements on the fringe 
patterns at two angles of incidence will be sufficient to enable the true nature of 
any close spectral complex to be deduced. The method will be applicable through- 
out the range of the instrument which will need no additional constructional 
accuracy. ‘The disadvantage of the method is that the auxiliary dispersive 
apparatus would need to be rearranged. 
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+ ABSTRACT. When making calculations on a circuit, containing both electric impedances 
* and transformers, it is frequently desirable to consider the transformers removed and the 
4 constraints they impose replaced by a rearrangement of the impedances connected to their 
§ terminals. Such ‘‘ equivalent circuits’? may not always be found ; the rules are here 
2 established for their formation, and also for checking, by inspection, whether the trans- 


9 former constraints are removable in this way, in any particular case. It is shown that the 


- equivalent electric circuit of a transformer, having any arrangement of magnetic paths, 
_ is derivable from its magnetic circuit by application of the topological principle of duality. 
‘This cannot be done if the magnetic circuit is non-planar, as in the case of a transformer 
possessing four or more windings with leakage couplings; a physically realizable circuit 
‘does not then exist. 

Under certain conditions the principle may be applied in reverse and the impedances in a 
| given electric circuit may be coupled together by a suitable transformer, so that the various 


% current and voltage constraints are unaltered. 


§1. INTRODUCTION 


HE simple geometrical relation which exists between inverse, or dual, 

circuits (Russell 1904) of electric impedances, based on the interchange of 

junction-points and meshes, is well known, and has been shown by Cauer 
(1934) to be based on the topological principle of duality. It is intended in this 
paper to show that an identical relation of pattern exists between the magnetic 
circuit of a transformer and its “equivalent electric circuit”. Furthermore, there 
exists a simple inverse relationship between the magnetic impedances of the 
transformer and the electric impedances in the equivalent circuit, in a manner 
analogous to the impedance relations inthe inverse electric circuits referred to above. 
By way of illustration, figure 1 shows such a pair of electric inverse (or dual) circuits. 
Given the first circuit A, a point is marked inside every mesh a, b, .... and one 
external to the circuit k; these points, when joined together as shown by the dotted 
branch lines, form the junction points of the dual circuit pattern B. ‘This dual 


circuit B now has impedances inserted in its branches Z,, Z,, Zs, .... which are 
proportional to the admittances aj, a, a3, .. . .in corresponding (crossing) branches 
ot the first circuit A, in the sense Z, : 2,72, :....=@, id) :@3:.... oF 


Z,,= K*a,, where K is any real constant. ‘This network transformation replaces 
mesh connections in one circuit by star-point connections in the other and 
vice versa. The element kinds must also be changed, thus : capacity C in Ist 
_ circuit is replaced by inductance L in 2nd circuit and vice versa ; resistance R is 
replaced by conductance G and vice versa. Corresponding impedance magnitudes 
in such a pair of circuits are thus reciprocally related. ‘This type of dual relation- 
ship is reversible, and starting with a circuit such as B the dual circuit A would be 
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derivable. Using the symbol == to mean “‘is dual to”’, the relationships be- 
tween such a pair of circuits may be summarized thus : 


La C- «R= = Coeds en eee (1) 


Again it will be found that voltages across branches in one circuit vj, Vg, .... are 
proportional to the currents in corresponding (crossing) branches of the other, 
dual, circuit, so that 

SST OO ee ae (2) 


If one of the circuits forms a non-planar figure (so that it cannot be drawn ona 
flat sheet of paper without any one branch crossing another) there exists no dual 
circuit, except by the use of a special artifice given by Bloch (1946). 

Analogous relationships will now be established for transformers and their 
equivalent circuits, taking for granted the rules for the setting up of dual mesh and 
branch patterns illustrated by Figure 1, and without employing the specialized 
notation and language of topology. 


PDI DIO DIMOND A Siz 
2 2 & & 4, of k 
aa ene 
h on 


naan 
Circuit A. Circuit B. (a) (0) 

Figure 1, Method of constructing dual circuits. Figure 2. Ideal transformers and the 

apparent connection of their 


windings into “equivalent cir- 
cuits ”’ 


(a) mesh-type transtormer; 
(b) junction-type transformer. 


§2, EQUIVALENT CIRCUITS OF IDEAL TRANSFORMERS -: 

Figure 2(a) shows an ideal 4-limb transformer having windings with an 
arbitrary equal number of turns, N, wound on a core of infinite permeability. The 
limbs, which may have different cross-sections, carry fluxes ®,, ®,, .... which,. 
meeting at a junction, must add up to zero; hence ©,+0,+0,+@,....=0. 
Differentiating and multiplying by N shows that the induced ee aon also. 
add up to zero : 

65 Peg Pe, Hehe = OF ee eee (3) 


Furthermore, the M.M.¥.s of the various limbs are in parallel and must be equal, 
since the transformer is ideal : . 


ijmly =i, es sade ere, ehenslis (4) 


These principles apply to ideal transformers having any number of parallel limbs 
carrying coils of equal turns. 

From equations (3) and (4) the transformer coils may be considered to be 
connected in an equivalent circuit as shown in Figure 2(a); the various coil terminal 
voltages are in series and a common mesh current 7 flows. Such an arrangement 
is not a true “equivalent circuit” in that the impedances of the various coils are 


i 


| 


y ( 


14 
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‘unknown since the original transformer was ideal. This circuit should, at this 
i stage, be regarded merely as a possible pattern into which the various voltages and 
currents of the transformer windings may be arranged so that the same constraints 


between them apply as are represented by equations (3) and (4). 
A second type of ideal transformer, shown in Figure 2 (6), has a series magnetic 
circuit; the equations for its currents and voltages are 


tigtizty+....=0 Se eee) 


-and Cp bay lg eC ay AL bah onente (6) 


so that the coils may again be imagined to be connected in an ‘“‘equivalent circuit”’, 


as shown in the lower part of Figure 2(b). Here the terminal voltages are all in 
» parallel and the currents entering the junction point must total zero.* © 


It may be seen that these circuits can be formed by applying the rules for dual 
figures, explained in $1, so that meshes in the first circuit are replaced by junction 
points in the second. In this present ideal case the equivalent electric circuit is 
the dual figure of the magnetic circuit, as has been illustrated by the superposed 
dotted circuit-in Figure 2 (6). 

The physically dual quantities are again those which lie on intersecting 


' branches but bear no direct relation to the dual electric impedances of equations (1) 


and (2). Thus, flux rate of change is dual, here, to voltage, while current is dual to 


M.M.F. in the intersected limb.t 
There again exists a constant of proportion between the quantities in such dual 


circuits. The duality of these electric and magnetic circuits is clearly reversible. 
Thus if we had started with the electric circuit we could have obtained the trans- 


former magnetic circuit by drawing the dual geometric figures (shown dotted in 
Figure 2 (b)) and replacing voltage e by flux rate of change d®/dt, and current z by 
M.M.F. M. ‘Thus the duality may be expressed as 


e=—=d0/dt im MS (7) 


One further example is needed to illustrate the complete rules for ideal trans- 
formers, viz. a transformer containing both series and shunt limbs. The trans- 
former shown in Figure 3 (a) is of such a type and, to make the case more general, 
two windings are shown on one of the limbs. All coils have N turns. With 


the fluxes as shown, the relation holds : 
OD, = (O;4+0,)=—(@3+G). tn eee (8) 


The schematic magnetic circuit of M.M.F.s (M,M,....) and fluxes has been drawn 
in Figure 3 (b), which also shows the dual figure (dotted lines) drawn according to 
the rules of interchange of meshes and junction points. On the assumption that 
this figure gives correctly the ‘‘equivalent circuit” of the transformer, it has been 
redrawn in Figure 4 with the dual currents and voltages marked in. Here the new 
junction points a, b....k have been identified. 


* The directions of the currents in the transformer windings and in the equivalent circuit may be 
related, if the sense of the coil winding is taken into account. Thus in Figure 2, with the windings 
as shown, the currents in them are assumed positive in the directions marked ; in this case the 
currents in the equivalent circuits have been made consistent. It is similarly possible to relate the 
directions of flux rates-of-change and of voltages. ne 

+An arbitrary set of rules could be applied for relating the directions of the currents and of the 
M.M.F.s in these two superimposed figures. For example, with reference to Figure 3, it could be 
stated that current direction is positive when it crosses a magnetic branch containing an M.M.F. 


which is positive on the left-hand side. 
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That this circuit is correct as regards the two coils on the limb carrying flux ®, 
follows from the fact that the voltages of the coils are identical, so that they may be 
regarded as being in parallel with one another. Again, the two limbs carrying 
fluxes ©, and ®, are magnetically in parallel and their M.M.F.s must be equal, or 
M,=M,,; the dual electric branches, a, b and b, k, carry equal currents, being in 
series. 

However, the branch a, k must be shown to be equivalent to that dual limb on 
the transformer which carries the flux ®,. The voltage across this coil is Nd®;/dt, 
which, from equation (8) is 


es = Nd®,/dt = N(d®, /dt + d®,/dt) = — N(d®,/dt + d@,/dt) =(e, + e2) = —(e3 +4), 
Ppl (9) 


and this is seen to be consistent with the mesh arrangement of the equivalent 
circuit in Figure 4. 
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(a) (0) (a) 1) 
Figure 3. The dual of a multi-limb Figure 4. (a) The ‘‘ equivalent circuit ”’ 
__ transformer. of the ideal transformer in Figure 3; 
(a) multi-winding ideal transformer; (6) the current and voltage pattern. 


(6) magnetic circuit (full lines) and 
the dual figure. 


Other examples may be treated in the same way. No matter how many limbs 
the transformer may have, or what may be the arrangement of the coils, their 
connection into an equivalent circuit of the type described may be derived by this 
routine method of drawing dual figures, provided that : (1) all coils have equal 
turns, NV; (ii) the transformer magnetic circuit is a planar figure. Similar relations 
exist between these dual electric/magnetic quantities as between the ordinary 
electric circuit dual quantities (voltage and current) which have been represented 
by equation (2). The fluxes (and hence their rates of change) in the various limbs 
bear certain ratios to one another; the voltages in the (dual) equivalent electric 
circuit have the same ratios. That is, in the steady state, 


JOD erGads < JeDerty,s one Oy Lela Bae fare = ee eee (10) 
A similar relation between the M.M.F.s and currents may be written 
Mi 2Mg 2 Mises spain =y oh ta cite hare eee (11) 


The question of magnetic and electric impedance duality analogous to (1) does 
not arise in these ideal transformers, but will be dealt with in §4. 


$3) INCLUSION OF) DHE EXTERNAL CIRCULE 


When a transformer forms a component part of a circuit it exercises certain 
constraints on the currents and voltages. These constraints arise by virtue of the 
magnetic couplings between the transformer windings to which the external 


en 
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; | circuit is connected. When the transformer is replaced by an equivalent circuit, 
i) this must be arranged so as to exercise the same constraints. In the case of the 
> ideal transformer this equivalent circuit has been derived by applying rules of 
iw! duality to the transformer magnetic circuit; the external loads and generators 
++ must not be included when these rules are applied, but must be imagined dis- 
connected from the transformer windings and then, when the dual figure of the 
_ magnetic circuit has been formed, they may be reconnected to corresponding 


points on this dual circuit so as to support the same currents and voltages as 


originally. 


For example, Figure 5 (a) illustrates an ideal transformer connected to a set of 


| loads and a generator E. This particular transformer has already been used as an 


example (Figure 3) to illustrate the use of the duality rules, from which it has been 
‘shown that the various windings may be considered to form the equivalent circuit 
‘connections in Figure 4. This “‘circuit”’ is more correctly to be regarded as an 
-arrangement of the various windings, currents and voltages into a pattern which is 
‘consistent with the constraints exercised by the transformer. It may alternatively 
be illustrated by a pattern such as that in Figure 4(4). The various loads Z,, Z,,.... 
and the generator E may now be connected in this pattern so as to support 
‘the correct voltages and currents. It will be noticed that the resulting circuit 
(Figure 5 (6)) contains no elements related to the transformer itself since, in this 
-case, this is ideal and possesses no leakage inductance. 


(a) (0) 


Figure 5. A loaded ideal transformer and Figure 6. A transformer with inter- 
its equivalent circuit having iden- connected windings and the failure 
tical current and voltage relations. of the “‘ equivalent circuit ’”’. 


‘(a) loaded transformer; (+) equivalent 
circuit. 

It may sometimes be found that it is not possible to carry out this process as, 
for example, when there is an external impedance common to two windings. 
Thus Figure 6(a) shows a transformer with two secondary windings having equal, 
induced voltages, e,, and common connections between them. ‘lhe two secondary 
currents are 7, and7z,._ The equivalent voltage and current “‘circuit”’, or pattern, 
formed by the dual figure of the magnetic circuit, is given by Figure 6(6). In this 
pattern the two secondary voltages naturally appear in parallel and there is no way 
in which the external circuit Z,, Z, and Z,’ may be connected so as to support the 
-same currents and voltages as in the actual system. In general it is only when each 
winding is connected to a separate load or generator that connection of such loads 
in an equivalent circuit is guaranteed. The equivalent circuit pattern merely 
represents one possible self-consistent arrangement of the transformer winding 
‘currents and voltages; an entirely different set of constraints may apply to the 
external circuit. In any given example it is usually possible to check the existence 


of an equivalent circuit by forming the dual figure of the transformer magnetic 


circuit and attempting to reconnect the external electric circuit to this figure, 
paying due regard to directions of currents and M.M.F.s. 
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§4. THE PRACTICAL TRANSFORMER—ITS MAGNETIC CIRCUIT 


In practice; transformer cores possess finite permeability, and their windings. 
set up leakage fluxes. A practical, loaded, transformer may be regarded either as 
an electric circuit having a certain arrangement of voltages, currents and imped- 


ances of as" magnetic circuit of M.M.F.s, fluxes and reluctances. It is intended 
now to show that these circuits are duals, as was the case for ideal transformers. 

Figure’7(a) illustrates a simple two-winding transformer with primary’and 
secondary leakage fluxes ®/ and ®{ anda finite core reluctance S, which may 


possibly include effects of air gaps. ‘The schematic magnetic circuit is\given in 
Figure 7(b). - S, and S, are the reluctances of ee (air) paths of the San and 
secondary leakage fluxes. 

The current i, flowing in’ the Nath primary coil sets up an M.M.F. Ni, 
which should. be regarded as the driving M.M.F. M,. 'The induced current’, in 


the secondary winding gives rise to a response or sue -u.u.F., M,.- ‘The resultant 


of these M.M.F.s sets up the working c core flux QD. 
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k indicates reluctance. (a) ; (0) 3) 
(c) eee 
Figure 7. The practical transformer. Figure 8. Equivalent circuit with: 
Its magnetic circuit and (dual) lossy magnetic material. 


equivalent circuit. 

(a) transformer with leakage fluxes; 

(6) magnetic circuit ; (c) equivalent 

circuit. 

The dual figure of this magnetic circuit (shown dotted) may now be formed. 
according to the rules, but the question arises : what electric elements are dual to the 
reluctances? Reluctance causes a loss of M.M.F. proportional to the flux passing 
through it, 4 =@S; in the (dual) equivalent circuit M.M.F. is replaced by its dual, 
current, and this must flow through an electric admittance to give a drop in voltage, 
dual to d®/dt. ‘This admittance is then the dual of reluctance, and furthermore 
it must be inductive. . 

This may be seen more clearly by dividing corresponding sides of equations 
(10) and (11): : 
Mie 1 te. tg , 
jo®, . jo®, . jw®, seen Se =e 


En SEOs 
that is, the ratio between the various magnetic impedances is the same as the ratio. 


between the corresponding electrical admittances. If the electrical elements be 
inductances this last equation may be written - 


My : M, é Ms; . pa i : 1 . 1 
We jady sae 


jol, " jol, Joly 
so that the duality relation is independent of frequency ae 
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_ It is seen that once again there exists a reciprocal relationship between the dual 


physical elements, S,, and L,,, just as in the case of electrical duality, where the 
dual impedance elements Z,, and 1/a, were reciprocals (§1), since _ 


S,=Ni,/®, and—-L,=N®,/j,, or S,L,=N?. ......(12) 


The equivalent circuit for the practical transformer, in the example, has been 
drawn in Figure.7(c). .This equivalent circuit is, of course, well known. The 
values of the inductance elements are given directly from the values of the dual 
magnetic reluctances from which they were derived : L, is dual to S, and so 
o D/(z,—7,); L, is dual to S,, and so « @j/i,; Li is dual to S,, and so « O{/i,. 
Note that this equivalent circuit, obtained by the duality rules, is still correct if the 
secondary be open-circuited or short-circuited. In the first case 7,=0 and the 
M.M.F. JM, vanishes; that is, it is magnetically short-circuited. In the second case 
e,=0 and the flux (®{—®) vanishes; that is, ®{=@ and the magnetic branch 
(containing M,) is open-circuited. 

The equivalent electric circuit of any other practical transformer is related to its. 
magnetic circuit in this same dual way, whatever the arrangement of the limbs, 
provided that the magnetic circuit (including leakage flux paths) forms a planar 
figure. 


§5. PRACTICAL TRANSFORMER—INCLUSION OF CORE LOSS 


To be correct practically, the magnetic circuit of a transformer must include 

loss elements which are due to the core hysteresis and eddy currents. The values. 
of such core losses and effective permeabilities are known for many common types. 
of magnetic material (Macfadyen 1947). 
The presence of hysteresis loss in a magnetic sample is indicated by the finite 
area of the hysteresis loop; the average loop is of such a shape as to necessitate a 
non-linear relationship between 8 and §, but if the production of harmonics be 
ignored then, as far as fundamental components are concerned, the hysteresis. 
loop may be regarded asathin ellipse. This elliptical shape is actually approached 
in practice if the effects of eddy currents are appreciable. Thus both Sand §- 
may be taken as having sinusoidal waveforms slightly displaced in phase, a point 
of view which has been adopted by Macfadyen (1947) in the development of the 
idea of complex permeability. 

It is convenient to use this idea here and to regard the magnetic circuit elements, 
reluctance and ‘“‘magnetic loss”’, as complex magnetic impedances, SS. 

The simplest magnetic circuit that may be taken for illustration is a ring having” 
a mean magnetic path length / and a cross-section A (Figure 8(a)). The reluctance 
of this core, Ss is related to the permeability, 1, by the expression S=(//A)(1/). 
amp-turns/weber (henries~+) = (//A)(r +7). These real and imaginary series 
components of reluctance are illustrated in the magnetic circuit of figure 8 (6). 

If a coil of N turns be wound on this core, carrying a sine-wave current, z, and 
having an induced sine-wave voltage, e, the magnetic field-strength and flux~. 
density become 


Ce Ni i ieamp-turns/miy ln 4) ewes (13) 


B= —jelwNA weberjm?, ~  —...... (14) 


assuming that the %/ loop is elliptical. 
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Dividing corresponding sides of (13) and (14) relates the magnetic reluctivity 
to the electric admittance : 1/f@=(7+j)=jN?’wA/lZmetres/henry. Thus an 


equivalent electric circuit may be found for the coil, shown in Figure 8 (c), which. 


‘includes the effects of the complex reluctance S= (1/A)(r + jis) =j N2w(G —7B) 
-amp- turns/weber(henries~!). Here the electric circuit is given as an admittance 


1 Ze (G—jB) shown in the figure as a resistance Rin parallel with an inductance L. 
Comparison of the real and imaginary parts of the above equation for S shows 


that : (i) the real part of the reluctance cg gives rise to a susceptance B, in the 
equivalent circuit : /7/A « B=1/wL (as has been shown already for the loss-less 


core) ; (ii) the imaginary part of the reluctance ‘S gives rise to a conductance G, and 
(2b/A) x G(=1/R). 

Notice that this equivalent circuit depends on angular frequency w and so may 
be applied only to transformers working at a constant frequency. 

This electric circuit, Figure 8 (c), is the dual of the magnetic circuit, Figure 8 (4) ; 
the two electric elements in parallel correspond to the two magnetic elements in 
series. However complicated the magnetic circuit of a transformer may be, such 
an equivalent electric circuit for each limb may be derived. 


§6. THE MULTI-WINDING PRACTICAL TRANSFORMER 


It is interesting to apply the method of duality to setting up the equivalent 
circuits of multi-winding transformers and to show that, in the general case, there 
can be no physical circuit if the transformer has four or more windings. 

Figure 9 (a) shows a transformer with one primary and two secondary windings, 
assumed to be coupled not only by the main working flux ® but also each to each by 
leakage fluxes. A possible fourth winding has been shown by dotted lines. In (4) 
the equivalent magnetic circuit is shown; here the M.M.F.s M,, M,, Mg, .... are 
those due to the currents in the primary and in the secondary windings. ‘These 
various M.M.F.s are shown coupled by leakage reluctance paths which are denoted 
by suffixes : thus S, is a reluctance shunting M, whereas S,,. couples M, and M, 
andsoon. ‘lhe addition of the fourth winding to the magnetic circuit is shown by 
the dotted lines ; a new secondary M.m.F. M; is coupled to the others by reluctances 
apg eae te: 

Considering the 2-secondary case only, the equivalent electric circuit may be 
-constructed dually in the usual way, by marking in reference points in every 
magnetic mesh and joining up to form the dual circuit (Figure 9(c)). The various 
inductances in this circuit are marked by suffixes corresponding to the various 
appropriate (dual) leakage reluctances, their respective magnitudes being given 
by equation (12). 

In the case of three secondaries, however, such a dual equivalent electric 
circuit cannot be formed because the magnetic circuit is anon-planar figure. The 
case considered here includes every possible arrangement of mutual couplings of 
the windings by leakage fluxes; equivalent circuits of multi-winding transformers 
have, however, been developed using negative elements. (For example see 
Blume, Camilli, Boyajian and Montsinger, 1938.) That this figure must be 
non-planar follows from the topological fact that five or more points in a plane 
(i.e. the M.M.F. “terminals’’) cannot be joined each to each by lines (i.e. leakage 
reluctance paths) without at least one cross-over. Such an essential cross-over is 
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si) marked * in Figure 9(b). This fact depends only on the number of M.M.E. 


“terminals” and not on how they are arranged in the magnetic circuit. 


| §7- TRANSFORMERS WITH UNEQUAL TURNS ON THEIR WINDINGS 


If the principle of duality is to be applied to the construction of equivalent 
| circuits for practical transformers, those cases cannot be excluded in which the 
_ turns ratios between the various windings are not unity. 

The duality relationship which exists between the magnetic and the electric 
| circuits of a transformer has been established in the preceding sections by starting 
with the idea] transformers of the “mesh type” and “junction type”’, illustrated. 


(a) 


© 
Figure 9. (a) A transformer with three Figure 10. (a) A bridge circuit and 
windings ; (b) its magnetic circuit ; (dotted) its dual figure; (b) the 
(c) the equivalent electric circuit dual transformer and its loads. 


and (dotted) the effect of adding a 
fourth winding. 


in Figure 2; the equations (3), (4), (5) and (6), which give the current-voltage or- 
M.M.F.—flux relations, involve the necessary condition that all the coils have an. 
equal number of turns, N. If the coils have an unequal number of turns, these 
equations are modified. ‘Thus if, in the mesh-type transformer, the limb fluxes,. 


®,, D,, D,, .... link with coils having, respectively, N,, No, Nz, .... turns, then 
equations (4) and (5) modify to Nyz,= Not2= N3i3= ...., that is all limb M.M.F.S,. 
being in parallel, are equal, and similarly e,/N,+ €,/N2+e3/N3+ ....=0, since the: 


sum total of the limb fluxes must be zero. 
The above equations may be divided by n, the highest common factor of N,, 


Nz, Nz, .-.-, giving 


n n n N,/n N,/[n  N,/n 
These equations represent the voltage and current relations that exist in various. 
“sections”? of the individual coils, where one ‘‘section”’ is assumed to have n° 
turns. Then, treating each section as a separate coil, the dual figure may be 
constructed and a nominally equivalent circuit to the transformer obtained. Such 
a circuit cannot, however, be the true equivalent circuit, because the n-turn. 
sections into which the coil on a particular limb has been divided are electrically in. 
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series on the limb but are in parallel in this nominally equivalent circuit; the 


voltages across the transformer windings, which are proportional] to the turns 
N,, No, Nz, ..--, do not appear in this equivalent circuit at all. 

Again, a junction-type ideal transformer may be treated in a similar way, and 
in this case, too, no truly equivalent circuit can be derived. 

In the case of practical transformers, with cores of finite permeability (possibly 
complex) and having leakage fluxes, true equivalent circuits consisting entirely of 
practical electric elements cannot be formed unless all windings have equal 
numbers of turns. However, by the inclusion of zdeal transformers in these 
equivalent circuits the usefulness of the duality relationship may be extended to 
transformers with non-unity turns ratios. Thus if such a practical transformer 
has windings possessing unequal numbers of turns N,, Nz, Nz .... these may be 
reduced, in a theoretical way, to an equal number u by associating an ideal two- 
winding transformer with each separate winding ; these various ideal transformers 
must then have turns ratios of N,/n, N2/n, . 

The equivalent electric circuit may now be found by the duality rules in the 
manner described in the earlier sections, though care must be taken to include the 
ideal transformers as part of the external loads and not to apply the duality rules to 
them. 


§8. INVERSION OF THE MAGNETIC-ELECTRIC DUALITY 
PRINCIPLE 


The duality principle has so far been applied to the production of an electric 
circuit equivalent to a transformer (i.e. a magnetic circuit) connected up to external 
loads and sources of power. As with other dual relationships, such, for example, 
as the electric circuit duality described in § 1, the procedure may be inverted and, 
given a planar circuit of purely electric impedances and sources of power, an 
equivalent transformer may be constructed. ‘That is to say, these electrical 
impedance elements may be assembled on to a transformer in such a way that the 
constraints between the various currents and voltages are unchanged. ‘That such 
an inversion of the process is possible has been shown in §2 and expressed by 
equation (7). 

By way of example, Figure 10 (a) shows a bridge circuit of any linear complex 
ympedances Z,,'Z5, Zs,,..-. Ze, Catryme branch currents ?,, 2), 25, .2-- tn ene 
source of power is included, for generality, in series with Z,. The dual figure 
may be formed by marking reference points inside every mesh a, b, c, .... with 
one outside, k, and joining up by lines (shown dotted) intersecting every branch of 
the circuit. ‘This dual figure gives the required arrangement (Figure 10 (d)) of the 
transformer limbs and windings which are loaded by the same impedances, 


Z,, Z,,.... Zs. In this transformer the various limb fluxes are dual to the circuit: 


branch voltages (more correctly, the flux rates of change if non-steady-state is 
considered) while M.M.F.s are dual to the branch currents in corresponding 
(crossing) branches. 

Such a transformer is in general] ideal, the core having infinite permeability 
and the number of turns on every winding being arbitrary, but equal. A practical 
transformer may only be derived 7f the original circuit contains inductances so that 
the various limbs, dual to the circuit branches, may possess finite reluctances. 

Obviously, if the original electric circuit is a non-planar one, possessing 
cross-over branches, it has no dual figure and, therefore, cannot be assembled on a 
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transformer in this way. However, part of the circuit may be planar, and those 


particular branches may be dealt with in the manner described and replaced by an 


equivalent transformer. 
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ABSTRACT. An apparatus for the investigation of the transmission of sound along 
filaments at frequencies between 1,000 and 6,000 cycles per second is described, both for 
unstrained specimens and whilst they are being elongated at a constant rate of increase 
of strain. Measurements of the dynamic elasticity and damping factors of filaments of 
polythene, neoprene and nylon have been obtained, and the correlation of these results 
with the molecular re-arrangements which take place during stretching is discussed. 


§1. INTRODUCTION 

HE properties of long chain polymers differ in many ways from materials 

composed of. smaller molecular units and, in particular, the mechanical 

behaviour is considerably modified when the molecular chains become 
oriented. Thus the stress-strain curve for such materials is complicated by 
orientation produced on applying stresses. In general the elastic moduli are 
found to be much greater for highly oriented materials than for such materials 
in the isotropic state. 

Also, with most of these substances if the stress is applied slowly, irreversible 
flow takes place, and the strain is a function not only of stress but of time.. 

The ordinary stress-strain curves of such materials do not, therefore, easily 
yield information about the molecular processes which take place during 
deformation. 

Some workers have applied alternating stresses of various frequencies to 
high polymer specimens and measured the resulting strains in efforts to obtain 
more complete knowledge of the mechanical behaviour of these materials. 
Recently, measurements of sonic and ultrasonic velocities in fibres of various long 
chain materials have been made in the U.S.A. by Ballou and Silverman (1944) 
and Nolle (1947). The work described here extends the scope of this type of 
investigation by determining the behaviour of filaments under conditions of 
steadily increasing strain. It was hoped in this way to separate the “elastic” 
and “plastic”? components of strain, and so help to correlate the mechanical 
properties with the molecular re-arrangements which are known to take place 
when these substances are deformed. 
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The method used was to send a sonic wave down a filament of the material,. 
and to measure the velocity of propagation, V. From this the dynamic value 
of Young’s modulus, £, could be determined, using the relation V=4/(E/p) 
where p is the density of the material. At the same time, the attenuation of the 
wave gives a measure of the “‘internal friction”’. 

In order to measure the variation of these constants with the amount of 
orientation present in the material, the filament was stretched at a constant rate 
of increase of strain, and observations made during extension. 

It may be shown (e.g. Rayleigh 1896 a) that for the propagation of longitudinal 
oscillations along a cylindrical specimen, if the wavelength of the oscillations. 
is large compared with the diameter of the cylinder, as in these experiments, 
then the wave propagation is essentially planar and the linear modulus, £, is. — 
applicable. ; 

Since the direction of propagation is along the direction of orientation, the 
anisotropy will not affect the nature of the wave motion, the relevant modulus. 
being that along the filament. 

§2. APPARATUS 

The apparatus (Figure 1) was designed to measure the velocity of longitudinal 
oscillations in filaments (about 1 mm. diameter) at frequencies between 1,000: 
and 6,000 cycles per second; a 
measure of the absorption of rf per =f) 
sound along the filaments could Rocoooocc 
also be obtained. 

The source of sinusoidal eee pins 
oscillations was a_ resistance 
capacitance bridge oscillator as 
described by Clifford (1945). 
The advantages of this circuit 
are great stability and high volt- 
age output. Variation of fre- 
quency was obtained by altering 
the capacities of a Wien bridge. 
These were ganged air conden- 
sers, and adjustment was made 
through a 6:1 reduction drive 
so that fine control was possible. 

5 ; Pre-amplifier 

The output from this oscill- eee 
ator was fed into a two-stage R.C. Output Double - beam connected to winding-drum 
coupled amplifier. The first ed 
stage was a Schmitt phase split- 
ter, and the second and output 
stage consisted of four 6L6 valves wired ina parallel push-pull. These valves were 
worked as triodes by connecting the screens through 100 ohm stoppers to the 
anodes and each valve was biased so as to work in Class A. The distortion 
produced was thus extremely small. A matched transformer coupled the 
output from this amplifier to the moving coil of the mains-energized speaker. 
As a check, both the current and voltage wave-form in this coil were examined. 
on a cathode-ray oscillograph, and no distortion from the sine wave input could 
be detected. 


M.C. Speaker &. 
Driving Cone 


~<— Filament. 


Crystal Pick-up | 


Figure 1. 
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Figure 2 (a) shows the construction of the Perspex cone which was an integral 
part of the coil former and which served to transmit the oscillations to the filament. 
_A stiff brass diaphragm carried this system in the centre of the speaker field. 

| By careful centring it was possible to ensure that only longitudinal oscillations. 
9) were produced. 

___ A crystal gramophone pick-up was mounted on an optical bench, fixed in 
| the direction of the filament so that the position of the pick-up could be varied 
j over a distance of about 14m. from the source of oscillation. A difficulty 
} encountered in early experiments was the great change in the cross-section of 
the filaments when extended to 200 or 300%, so that a fixed contact needle became 
; too loose for large extensions. To overcome this, a Perspex needle with a spring- 
* loaded contact was designed, whereby contact closed on the filament as the 
| cross-section decreased (Figure 2(6)). The small voltage output from the 


©) ~— Spring Balance 


Cord attached to filament 
SA 


Speech Coil 


Brass Diaphragrn 


~— Perspex Locking Nut 


(a 
Hole through which 
filaments are threaded 
Driving Shaft 
a=— Screwed Plunger AS Wading 
ie eae 
on Loading Spring 
7 ae a 
‘is Ball Bearings if Filament: Cord tension shown thus —>— 
Figure 2. (a) Perspex coil former and Figure 3. 


driving cone. 

(6) Perspex pick-up needle. 
| crystal was amplified by a single stage pentode amplifier located near to the 
} actual pick-up. This was of conventional r.c. type. The output from this. 
_ was then carried by a screened lead to a Cossor double beam oscilloscope, model. 

339, and fed through the two stage amplifier incorporated in this instrument 
* on to the Y, plates. The screening of the first amplifier was very thorough, as 
was the decoupling of it from the oscillator and amplifier. Even so, a smalf 
amount of 50 c/s. mains ripple was found to be present, and a three element 
high pass r.c. ‘T’ filter was placed in the lead from the pre-amplifier to the 
cathode-ray oscilloscope; this was found to be very effective in eliminating these 
lower frequencies. 

The voltage applied to the speaker coil was carried by screened lead to the 
other Y plates of the oscilloscope, no amplification being necessary. By adjust- 
ment of the time-base and synchronizing control, the trace length could be 
adjusted so that a single wave corresponded to any set distance on the tube. 
Phase difference between the two traces could then be measured as the distance 
between corresponding parts on the two sine waves. Any residual differences 
of phase, due to amplifiers, etc., were eliminated by the system of measurement, 
whereby a series of measurements of phase shift against distance along the 
filament was made. 

PROC. PHYS. SOC. LXII, 2—B 
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The free end of the filament was joined to a fine cord which was attached to 
the drum and tension measuring instrument at the far end of the bench (5 m, 
from the source). The principle of this instrument is shown in Figure 3. The 
spring balance was calibrated by loading the cord with known weights. The 
drum was driven by a }h.p. induction motor through a double worm reduction 
gear which could be varied, the system being effectively synchronous. The 
gearing used in these experiments was such that the end of the filament was 
drawn at the rate of 0-194 cm/sec. Whilst the filament was being stretched, 
photographs of the cathode-ray oscilloscope were taken at regular time intervals 
(corresponding to equal intervals of strain) with a 16 mm. ciné-camera, and the 
film measured later. This enabled much more accurate measurements of 
phase shift to be made than by direct observation. 

The frequency of the oscillator was checked against the supply mains. The 
oscillator was set on to a harmonic of the supply frequency, and by using the 
double beam oscilloscope, the exact harmonic was readily found and hence 
the frequency determined to the accuracy of that of the supply. 

Variations in the latter may have caused errors in the frequency of 1-2%, 
but these variations were of long period and, whilst causing errors in the absolute 
values, did not affect the relative results for any individual experiment. 

The voltage applied to the Y, plates of the cathode-ray oscilloscope was 
measured by a suitable rectifier-type a.c. voltmeter so that an estimate of the 
damping in the filament could be made. 


$3: LHEORY OF DHE MIELEoOD 

The system considered is a filament attached to a stiff massive diaphragm 
with the pick-up acting as a partially reflecting boundary at a distance / along the 
filament. The total length of the filament is considered sufficiently great for 
reflections from its other end to be neglected. It is assumed that varying the 
position of the pick-up will not affect the amplitude of vibration of the diaphragm, 
and measurements with the apparatus employed showed that this was so. 

Let the displacement of the diaphragm &, be given by 


éj= Asin @i.) Ya age «lS See (1) 
and the displacement of the pick-up by 
€ = Bsin(@t'G) "oe ee (2) 


where w = 27 (frequency). We now need to find B/A, the ratio of the amplitudes, 
and 6, the phase difference, in terms of the position of the pick-up, /, the reflection 
coeficient at the pick-up, m, and the propagation constants k and «, where 
k=w/V (V =velocity of sound) and « is the damping coefficient. 

During the initial transient stages, reflections of the progressive waves will 
occur, until a state of equilibrium is reached in which a system of stationary 
waves is set up in the filament. 

In the steady state the system can be represented by the sum of a single 
progressive wave travelling down the filament and a single reflected wave travelling 
back. 

Let the displacement of the outgoing wave €, at a distance x along the filament 
be given by &=ae~* sin(wt—kx); then the displacement €, of the reflected 
wave is given by = —mae°®- sin (wt —k(21—x)) and the resultant displace- 
ment € is given by 


E=£,+£,=ae* sin (wt — kx) — mae? sin (wt —Rk(21—«x)). 


_ At the pick-up 


: At the diaphragm 
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é, = ae" = m) sin (wt -- kl). ais veleiere (3) 


€,=a{sin wt — me sin (wt — 2k/)} 


=ay/{1 —2me-> cos 2k1 + me sin (wt —kI+B) «22... (4) 


where 


tan B={(1 me /(T—me) tan Ri ss wet (5) 


) Comparing (3) and (4) with (1) and (2) we see that 


B/A=(1—m)e™/4/{1 — 2me-™ cos 2kl + MiG} Sees (6) 
and 6 = given by (5). 


(a) Estimation of V, the velocity of propagation. 


The angle @ represents the phase difference between the speaker diaphragm 
and the pick-up and is measured as a function of /. Now 


tan 6 ={(1 + me) /(1—me-*™)tankl. ss a (7) 


For large values of «/ this can be written approximately tan 6 =tan kl whence 0 =i. 
The actual curve obtained from (7) is of the form of damped oscillations 
about the line 0=Al, cutting the line at points of inflexion where k/=0, 7/2, =, 


37/2 etc., the interval in / being A/4. Figure 4 gives an experimental curve for 


polythene showing how the curve approximates to a straight line for large values 
of]. Itcan be seen that an accurate determination of the gradient can be obtained. 
This gradient gives a value of k. Hence, since V=27f/k, the velocity may be 
determined. 


(6) Estimation of «. 
The expression for the amplitude B/A (eqn. (6) above) is more complex. 
A graph of B/A against / takes the form of damped oscillations (cf. Figure 5). 
‘The rate of decay and the number of oscillations present depended on both « 
and m. For large values of « such as occurred with neoprene, the expression 
can be simplified to B/A =e. This enables a direct value of « to be obtained 
from a logarithmic plot of amplitude against distance. 
For smaller values of « the following method was used. Equation (6) can be 
written 
Agee Co eine COS 2k1 
BI lane 


ol 2 


Z 
~ (l—m)? (1—m)? (1—m)? 


For large values of al the second term vanishes and the expression for A2/B? 
is of the form of small oscillations about an exponential rise. From (8), neglecting 
the second term for high values of «/, 


In (A?/B?) =201+C. 
‘Thus « is determined by the slope of the graph of In (A?/B?) against /. 
None of these methods could be employed during the drawing of the filaments 


as the data available from one measurement were insufficient to give a value of «. 
S=2 
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§4. EXPERIMENTAL PROCEDURE ? 

In carrying out an experiment the following procedure was adopted. One 
end of the filament was attached to the cone and the other end to the cord; the ' 
motor was then run until the filament was just taut. Measurements of phase ' 
and amplitude were then made with the pick-up, at centimetre intervals along — 
the filament. Figures 4 and 5 show the type of results obtained. From these 
measurements, both V and « for the unstretched filament could be determined. 


1440.— =] Tol 


1260 +- t If | 


Phase (degrees) 


Amplitude (arbitrary units) 


20 


0 20 40 60 80 


Distance (cm) 


eee 
0 20 40 60 80 100 0 


100 
Distance (cm) 


Figure 4. Polythene 3,000 c/s. Figure 5. Polythene 1,500 c/s. 


In order to measure the variation of the velocity of sound with strain, the 
pick-up was then fixed at a sufficient distance along the filament for the approxi- 
mate relation =k] to hold. ‘The motor was started and readings of the spring 
balance and photographs of the cathode-ray oscillograph trace taken at convenient 
intervals of time. 

In the case of nylon, where the damping was extremely low, the approximate 


form of the relation between phase and distance along the filament was found 
to be inaccurate and the accurate form . 


tan 0 = {(1 — me?) /(1 + me) tan kl = P tan ki 
had to be used. 

In general, values of « were not obtained whilst the filament was being 
stretched. In the case of neoprene, however, a separate set of experiments. 
was made in which the filament was stretched by various amounts, held there, 
and a separate determination of V and « made as above. This was possible 
for neoprene as the plastic flow was not so rapid as to invalidate the results. 


All the experiments were carried out at room temperature; this is stated with 
the results. 


§5. EXPERIMENTAL RESULTS 


(i) Polythene (polymerized ethylene). Results were obtained for filaments 
of polythene of approximately 0-6 mm. diameter. Prior to testing, batches of 
three lengths were placed in a glass tube five feet long, and kept in a long oven 
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j at a temperature of 90° c. for several days, and then allowed to cool slowly, free 
| from restraint, in a horizontal position. 

' The filaments were stretched at a rate of 0-0015 per second. Results for the 
)) velocity of sound and for the dynamic modulus of the material are given in Table 1 
| and the dynamic modulus is plotted against strain in Figure 6. From the values 
| of the tension built up in the filament during drawing, values of the true ‘static’ 
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| stress were obtained and are also given in Table 1. Values of the gradient of the 

_ stress-strain curve, the nominal ‘‘modulus’’, are plotted in Figure 6. A direct 
comparison between the values for dynamic and static moduli is thus given. 

In certain experiments the winding machine was stopped when the strain was 
1-0, and the relaxation of stress was then measured against time. ‘This is given 
in Table 2; the change in dynamic modulus during the relaxation of stress was 
found to be extremely small. 


Table 1. Polythene. 
(Room temperature 20°c. Frequency 3,000 c/s.) 


Stress Velocity pees Stress Velocity meee 
“Strain (dynes/cm? of sound ( Stale Strain (dynes/cm? of sound (dynes soma’ 

x 107) (m/sec.) sie 10%) x107) — (m/sec.) 10°) 
0-000 0-459 782 5-65 0-187 10-78 780 5-60 
0-014 2-25 782 5°65 0-227 Iilesy 786 70) 
0-020 2:60 784 SIH 0-267 11-95 797 5-82 
0-027 3°58 787 5-70 0-441 forts 863 6-88 
0-033 4-01 787 5-70 0-641 19-71 987 8-98 
0-040 4-82 - 782 5-65 0-774 22°85 1070 10°55 
0-053 5°72 780 5-60 0-908 25°8 1122 11-60 
0-067 6°53 777 5:56 1-04 28-78 1207 13-41 
0-081 7:24 770 5:47 1:11 30°55 1254 14-50 
0-094 7°81 770 5:47 1:24 33-50 1286 15-80 
0-107 8:51 767 5-44 i237 36°42 1320 17-20 


0-146 9°81 767 5°44 
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Most of the experiments were carried out at 3,000c/s. Some work with 
undrawn filaments was also done at 1,500 and 6,000 c/s. Values of the damping _ 
factor and velocities for these frequencies are given in Table 3. 


Table 2. Polythene. Relaxation of Stress 
(Strain 1:00. Temperature 20° c.) 


Time Stress Time Stress Aeane Stress 
2 2 2 
(seconds) oe (seconds) es (seconds) Ree: 
0 PED} 280 24-8 2560 22°55 
15 26°8 580 pier 4360 22-22 
40 26:2 760 23-58 8000 21-50 
65 25:8 1360 PPEA iT 64000 20-20 
90 25'S 1960 22-68 


Table 3. Polythene. Variations with Frequency 


Frequency (c/s) 1500 3000 6000 
Damping factor a (cm?) 0:006 0:012 0-033 
Velocity of sound (m/sec.) 745 782 811 


Table 4. Neoprene 
(Room temperature 20°c. Frequency 3,000 c/s.) 


Stress Velocity Dynamics Stress Velocity Dynan 
: modulus ; modulus 

Strain (dynes/cm? of sound (Hexeslon! Strain (dynes/cm? of sound Ca 

1103) (m/sec.) 10°) < 1102) (m/sec.) 108) 
0-10 0-50 133 2°359 3-36 27-61 790 82:7 
1-76 799 311 12-88 352 30:10 845 95-0 
1-92 9-15 352 16:5 SOF) 34:17 913 110°5 
2-08 10-36 412 DIDS 3°84 38-10 1000 133-0 
2:24 12-20 443 26:35 3°99 41-60 1017 137-5 
2°395 13-90 480 30°75 4-16 44-30 1052 146-2 
2°56 15-79 545 39°5 4-32 49-30 1110 164-4 
Dee, 17-85 588 45-9 4-48 53-40 1218 196°8 
2:88 20-20 634 53-4 4-64 58-50 1280 218-3. 
3-04 22°65 692 63-7 4-80 a 1430 Pegs: 
3'-19 25-30 726 70:3 


Table 5. Neoprene Damping Coefficient («) 
(Frequency 3,000 c/s.) 
Strain 0-06 0:20 0:40 0-80 1:00 1:20 1-40 1:60 1:80 . 2-00 
a (cm~) 0:51 0:47 0-44 0-40 0-29 0:19 0-18 -0:12 0:07 00-6 


(ii) Neoprene. ‘The second material investigated was a synthetic rubber 
Neoprene GN. It was found convenient to strain this material at twice the rate 
used for the other materials; the rate of straining was thus 0-003 per second. 
Preliminary experiments showed that this did not affect the results appreciably. 
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» As with polythene, the gradient of the stress-strain curve was calculated. 
| Figure 7 gives a plot of the dynamic modulus and the gradient of the stress— 
_ strain curve against strain; numerical results are given in Table 4. 


‘Some results for the damping coefficient were obtained as previously described ; 


4 these are given in Table 5. 


(ii) Nylon. Some experiments were also carried out with nylon filaments 
of 0-7 mm. diameter. For large strains, sharp discontinuities (‘necks’) appeared ; 
these introduced large variations in the diameters of the filaments and rendered 
results inaccurate. However, approximate values were obtained and, as before, 
these are plotted together with the gradient of the stress-strain curve in Figure 8. 


_ The damping coefficient « for the undrawn filament was estimated from the 


amplitude measurements and was found to be about 0-003 cm}. 

Some experiments were also carried out with drawn nylon filament, 0:3 mm. 
in diameter. The dynamic modulus at 3,000 c/s. for this material was found 
to be 5-8 x 10!° dynes/cm?. This is greater than the highest value that could 
be obtained by drawing the unstrained filaments in the apparatus, showing that 
under these conditions rupture occurred before orientation was completed. 
The damping was extremely low and no reliable estimate could be made with 
the length of filament available. 


§6. SUMMARY OF RESULTS 


The results may be summarized as follows: 

1. In the unstretched specimens of all these materials, the dynamic value of 
Young’s modulus obtained from the measurements of the velocity of sound 
was found to be several times as great as the limiting gradient of the stress-strain 
diagram. 

2. The values of the dynamic moduli for both polythene and nylon, as the 
filaments were stretched, first showed a slight decrease followed by a rapid rise. 
In the case of neoprene, however, the dynamic modulus increased steadily. 

On the other hand, the gradients of the stress-strain curves for polythene 
and nylon decrease continually with strain, showing that the effects of plastic 
flow completely mask the changes in elastic behaviour. For neoprene, the 
gradient increases but at a much slower rate than the dynamic modulus. The 
latter had increased more than a hundredfold when the specimen had reached 
5-8 times its original length. | 


3. The results for polythene at 1,500, 3,000 and 6,000 c/s. showed that the 
damping increases with frequency, whilst there is comparatively little rise in 
the velocity of sound, and hence in the dynamic modulus. 

4. The damping coefficient of unstretched neoprene was found to be very 
much higher than those obtained for polythene and nylon. Up to a strain of 
0-8 the value of the damping coefficient showed little change; for strains greater 
than this, however, it decreased rapidly and had reached about one-eighth of 
its original value at a strain of 2-0. 


§7. CONCLUSIONS 


The above results conform, on the whole, to those expected from the molecular 
models already postulated for these materials from x-ray and other data. 
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Neoprene, like natural rubber, in the unstretched state, is generally assumed 
to consist of long molecular chains in a curled-up configuration; when a stress 
is applied the chains uncurl (Meyer et al. 1932, Guth and Mark 1934, Kuhn 1936). 
"The very low value of Young’s modulus under these conditions may be attributed 
to the ease with which rotation may take place around single bonds. When the 
chains become uncurled, further extension is only possible by changing the bond 
angles and this will correspond to a very much higher modulus. 

The rapid increase in the dynamic modulus with strain may thus be expected 
to correspond to a decreasing possibility of rotation around single bonds. ‘The 
large drop in the damping coefficient with extension, previously noted by Cassie, 
Jones and Naunton (1936) for slow oscillations, might also be expected from the 
more rigid molecular arrangement in the strained state. 

From x-ray data, Bunn (1941, 1942) has deduced that the structures of both 
polythene and nylon consist of long chains with small crystalline groupings, 
each of which is made up of portions of several molecular chains. Inthe unstrained 
material these have quite random orientations. As the material is stretched 
the chains become aligned and the crystalline groupings orient along the direction 
of stretching. When the alignment is completed the material resists further 
drawing and again behaves in an elastic manner for small extensions. 

With the polythene specimens available, it was found impossible to cover 
the complete stress-strain curve, as the filaments broke before cold drawing was 
completed. However, measurements were made with unstressed specimens 
of completely cold drawn nylon filament, which showed that the dynamic elastic 
modulus was about three times as great as for the undrawn filaments. 

The reason for the initial decrease in the value of the modulus, on stretching, 
in both these materials, is not altogether clear, but may be associated with the 
breaking up of random crystalline groupings on stretching. ‘Thus, for small 
strains, the molecular chains may be able to move more freely than in the unstrained 
material. Further increase in the strain will cause new groupings to be formed 
with their axes parallel to the direction of stretching. 

The frequency variation of the damping coefficient for polythene is in qualita- 
tive agreement with that expected from theory (Rayleigh 1896b, Bergmann 1938). 
The frequency range covered, however, was too small to derive any exact quantita- 
tive relations. At higher frequencies the period of the oscillations would 
become comparable with the relaxation times of the molecular chains, and it is 
intended to investigate this by extending the frequencies to these higher values. 

The large divergence between the static and dynamic moduli in the unstretched 
filaments may be due to similar effects and to the adiabatic conditions existing 
in the dynamic measurements. 

The thermodynamics of long chain systems has not been fully investigated, 
and it is hoped that a study of the dynamic properties of these materials under 
varying conditions of temperature and frequency will add to our knowledge of 
this subject. Work is now continuing on these lines. 
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A Theory of the Angle of Contact 


By G. D. YARNOLD Ano: B. J... MASON 
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_ ABSTRACT. It does not appear to be possible to give a quantitative account of the angle 
'-of contact of a continuous liquid in terms of cohesional forces. By using Young’s energy 
relation and making allowance for the effects of adsorbed films on the surface of the solid, 
-expressions are derived for the advancing and receding angles, and a general explanation of 
the hysteresis of the angle of contact is put forward. 


§1. INTRODUCTION 


s is well known, there exists no satisfactory theory of the angle of contact 
A between a liquid and a solid. Qualitatively, the existence of an angle 
of contact has been accounted for in terms of cohesional forces. Thus, 
Newman and Searle (1933) state that .an element of liquid lying in the line of. 
-contact experiences two cohesional forces, due respectively to the solid and the 
liquid, whose resultant must be normal to the tangent plane to the free liquid 
surface. Unfortunately, as is shown below, this conception cannot be made to 
yield a quantitative account of the phenomenon. 
Young’s equation expresses the cosine of the angle of contact # in terms 
-of the surface energies of the liquid—vapour interface, y;,y, the solid—vapour 
interface, y,, and the solid—liquid interface, y,. ‘Thus, 


CoSU=(Yi-Ye)ivry: ~<a (1) 
Recently, Doss and Rao (1938) have considered the effect on the surface energy 
-of the adsorbed film of liquid which condenses from the vapour so as partially 
to cover the free surface of the solid. | Assuming the fractional area o of the 
-solid which is covered by: the adsorbed vapour to be unaffected by the subsequent 
adyance of the bulk liquid over the solid, Doss and: Rao deduce expressions for 


-y, and y, and find... ... ears we 
cogf=2o—1, Meaawee? SAO. a6 ee ks (2) 
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Thus it appears from their theory that a zero angle of contact occurs only when 
the adsorbed layer completely covers the solid. Their concept of a partial 
adsorbed film, however, marks a definite advance in terms of which we have 
been able to give some account of the known dependence of the angle of contact 
on the state of the surface and the direction of motion of the bulk liquid. It — 
appears that the fractional area of the adsorbed liquid film is changed when the 
bulk liquid advances to cover the solid. 


§2. COHESIONAL. FORCES AND THE ANGLE OF CONTACT 

Following the classical conception of cohesional forces, consider a thin 
wedge of matter contained between semi-infinite planes inclined to each other 
at an infinitesimal angle 56. The attraction experienced by an element of mass. 
dm, situated on the edge of the wedge, due to the material composing the wedge, 
is directed into the wedge in a direction perpendicular to the edge, and may be 
written as kimSd, where k is a constant depending on the material composing 
the wedge. It follows that the force of attraction exerted by a wedge of homo- 
geneous isotropic material of finite angle « is directed along the bisector of the 
angle and is of magnitude 2kdm sina/2. 

Consider now the equilibrium of an element 5m lying at the point P in the 
line of contact of a liquid and a solid (Figure 1). The element 5m experiences. 
three forces Fy, Fi, and Fy due to the attraction of the solid, liquid and vapour 
phases respectively. Thus, 


F,=2k8m;  F,=2k,3msin0/2; Fy =2kySm cos 6/2, 


where kg, k;, ky denote the specific attractions of the solid, liquid and vapour — 
respectively and 6 is the angle of contact. ‘The condition that the resultant 
force shall be normal to the free surface of the liquid is then 

F, sin @ — Fy, cos 6/2 — Fy sin 6/2 =0, 
from which it follows that 2k, =k,,+ky. 

It appears, therefore, that the requirement that the resultant force acting ~ 
on an element dm should act normally to the free surface does not in any way 
determine the angle of contact of the liquid, but leads only to a highly artificial. 
result of no value whatever. We must conclude that no definite angle of contact 
could exist for continuous, homogeneous substances obeying a classical law 
of cohesion. ‘This result is in marked contrast with the comparable case of 
surface tension, of which a coherent account may be given along the lines first. 
laid down by Laplace. No doubt, the equilibrium existing in the line of contact 
between a solid and a liquid involves in reality the balancing of forces directed 
towards individual molecular centres. In the absence of detailed knowledge 
of such forces, however, we must be content with a more limited account of the: 
phenomenon in terms of Young’s energy equation. 


§3. SURFACE ENERGY AND THE ANGLE OF CONTACT 
Following Doss and Rao, we may suppose that the free surface of the solid. 
is partly covered by a liquid film, both before the bulk liquid advances to cover 
the solid and after the bulk liquid has receded. Moreover, if the solid surface 
is only partially covered by condensed vapour, then under normal conditions the- 
remainder is probably covered by an adsorbed film of air and other impurities, 
while under vacuum conditions a portion of the surface may be actually bare. 


Peewee 
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In general, therefore, the solid surface may be supposed to consist of a chequer 
pattern of small patches. As the liquid advances or recedes over the solid surface 


| its motion is easy over some patches and difficult over others, depending on the. 


nature of the actual molecular forces. The line of contact is generally irregular 
and the angle of contact variable from point to point. Indefinite spreading of 


_ the bulk liquid would not take place provided the patches of adsorbed liquid 


were in the nature of “islands” in the layer of adsorbed air. The angle of 
contact as measured by macroscopic methods is an average value whose cosine: 
can be expressed by Young’s equation, provided allowance is made for the 
fractional areas of the different patches composing the solid surface. In what. 


_ follows it is assumed that adsorbed films are of sufficient thickness for the interface: 


on either side of the film to possess the full interfacial surface energy characteristic. 
of the substance in bulk. 


Figure 1. Figure 2. Figure 3. 
Liquid advancing. Liquid receding. 

Let the surface energies of the interfaces solid—liquid, liquid—vapour, solid— 
air, liquid—air, air-vapour be denoted by yey, yiy, Ysa, Yta» Yav Tespectively. 
When the liquid is advancing slowly over the solid (Figure 2) let a fraction o, 
of the surface PA near to P carry condensed vapour, a fraction ¢, adsorbed air, 
etc., while a fraction 1—o,—¢, is bare. Similarly, let a fraction 4, of the com- 
posite solid—liquid interface PB near to P carry adsorbed air, while a fraction 
1—¢, is a true solid—liquid interface. We cannot assume that ¢.=¢), since 
the advance of the line of contact may tend either to compress or partially to remove 
the layer of adsorbed air. Once the bulk liquid has covered an adsorbed liquid 
film, no separating interface exists between them. ‘Thus, the average surface: 
energies of the interfaces PA and PB near to P may be written 


¥1= (Ys +Yxv) + ¢i(¥sa +Yav) +1 — 01 — da) ysv> 
¥2=$ol¥ga + Ya) + (1 — $2)y¥s3- 
For slow advance Young’s equation yields 


R Ys YAV y 
cos 0, =(o, +2 —1) aa oy eae + ($1 $2) =* +o, — ob, — +04, 
YLV YLV YLV YLV i 


where 6, is the advancing angle of contact. 
If it is supposed that the solid surface is perfectly degassed, then ¢,=¢, =0,. 


and equation (3) reduces to 


cos 0, =(1 —0;) (—) 0. eee ere ee (4) 
YLV 


124 G. D. Yarnold and B. F. Mason 


If, in addition, the solid is completely covered by a film of liquid condensed 
from the vapour in advance of the bulk liquid, o,=1 and the advancing angle of 
contact is zero. If, in the other extreme, no appreciable film condenses, o,=0 
and cos6,=(ygy—ysr)/yzy. This latter expression has been used in the 
past to define the so-called “equilibrium” angle of contact. According to the — 
foregoing theory, however, the expression defines an upper limit to the contact 
angle which corresponds to the advance of the liquid over a completely degassed 
surface and in the absence of any film of condensed vapour on the solid. The 
actual contact angle for a degassed surface must lie between these theoretical 
limits and is determined by the value of o,, and hence ultimately by the vapour 
pressure of the liquid. 

Similar expressions may be derived for the receding angle of contact, Og. 
For, let a fraction o3 of the surface PA near to P (Figure 3) remain covered by 
liquid and a fraction ¢, by adsorbed air, etc., when the bulk liquid recedes. Let 
a fraction ¢, of the interface PB near to P carry adsorbed air. Then the average 
surface energies of the composite interfaces PA and PB are : 


¥3=93(Ysr+ Viv) + bs(¥sa + Yay) + (1 — 93 — $3) ¥ sv, 
¥a=al¥ga + Yx4) + (1 — $a)ysu- 
For a slowly receding liquid Young’s equation gives 


cos 6; = (o,-+¢4—1) BE Hd asd Hbe—ty ng tae ay ee TO3% 
YLV MANY ny YLV YLV 


Again, if it is supposed that the solid surface is completely degassed, so that 
‘hs = b4=0, equation (5) reduces to 


cos 62 = (1 —as) Ga Ogee © Sl eee? (6) 
YLV 


If a complete film of liquid is left by the receding liquid, whether of macroscopic 
or microscopic thickness, o3=1 and the angle of contact is again zero. On the ~ 
other hand, if no film whatever remains when the bulk liquid recedes, then 
COS Op =(Ygy —ysz)/Y~tv- Thus the range of theoretically possible values of 
the angle of contact against a completely degassed surface is the same whether 
the liquid is advancing or receding. 


§4. HYSTERESIS OF THE CONTACT ANGLE 

Experiment shows that the advancing angle of contact between water and wax 
‘surfaces increases slowly with increase in velocity of the bulk liquid, while the © 
receding angle decreases both with increase in velocity and with increase in 
the time of immersion (Yarnold and Mason 1949). 

Under normal conditions, no part of the solid surface is likely to remain bare. 
Hence we may write o, +¢,=03+¢,=1. 

Let us now suppose that the advance of the liquid over the solid is accompanied 
by a partial removal of the adsorbed air film and a corresponding increase in the 
fractional area of contact of liquid and solid. Similarly, we may suppose that 
the recession of the liquid is accompanied by a decrease in the fractional area of 
contact of liquid and solid and a corresponding increase in the fractional area 
occupied by adsorbed air. Thus we may write ¢,-—¢,=« and ¢,—¢,=, where 
a and f are themselves functions of the velocity of the bulk liquid. During 
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j, immersion, we may suppose that there is a gradual replacement of air by liquid, 
the air going into solution. Thus ¢,—¢,=5, where 5 is a function of time. 
For the difference in the angles of contact equations (3) and (5) now yield 


Coe), cos0 =a — 6) (@ ma 4) =?) a5 (“2 asain =H) 
VEV. YLY 


If the time of immersion is very short, the possibility of the solution of air 
during immersion may be neglected. It follows then that ’ =0 and o, —o, =a —f. 
Thus equation (7) reduces to 


a—B 
(cos Og — C08 4), .9 = ae (Ysut+Yiv—Ysa—Yav)) oo yeeeee (8) 
Ly 


Provided y,y is small, it follows that an increase in the velocity of the bulk liquid 
_ is accompanied by a greater hysteresis of the angle of contact. At low velocity 
the receding angle appears to change more rapidly with change of velocity than 
does the advancing angle, suggesting that 8 varies more rapidly than « At 
high velocity it is known that the difference of the contact angles attains a constant 
value, but it is conceivable that no theory of this limiting case is possible without 
entering into hydrodynamical considerations. 
If the velocity of the bulk liquid is constant, the first term of equation (7) 
_ is constant, and the hysteresis of the angle of contact increases with increasing 
| time of immersion. 
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ABSTRACT. The effects (a) of the velocity of the liquid surface, and (4) of the time of 
immersion, on the angle of contact between clean surfaces of water and wax have been 
investigated by a dynamical method. It is clear that Ablett’s measurements are vitiated by 
the fact that the time of immersion is disregarded. The term “‘ equilibrium angle of contact” 
appears to be meaningless. 


§1. INTRODUCTION 
METHOD recently devised by one of us (Yarnold 1946) for the investigation 
aN of the hysteresis of the angle of contact has now been used to study the 
effects of the time of immersion and of the velocity of the liquid on the 
angle of contact between clean surfaces of water and paraffin wax. The range 
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of variation of the velocity of the liquid is necessarily quite small, extending 
only from 0-1 mm. per minute to 1-6 mm. per minute, but the results are of 
particular interest in that they discriminate clearly between the true velocity 
effect and the effect of the time of immersion in the case of the receding angle. 


§2. EXPERIMENTAL 
A plane water surface is raised or lowered hydraulically at a definite rate 
while it is in contact with the surface of a steel sphere coated with paraffin wax 
and suspended from a Sucksmith ring balance. Continuous observations are 
made of the force experienced by the sphere and of the position of the water surface. 
"The angle of contact is determined graphically (Figure 1) from the coincidence 


Times of Immersion 
1 3 mins. 
2 5 
3 75 » 
4 135 
5 38 


ADVANCING (a) 


Angle of Contact (degrees) 


0:10mm/min 


90 


| 
2 
3 
l | 3 
03 5 

-01 0 0-1 0:2 0 i0 20 30 40 50 

ft Time of Immersion (minutes) 
Figure 1. Figure 2. 


Receding velocity 1°6 mm/min. 


-of the observed upward force with the calculated hydrostatic upthrust on a segment 
of the sphere, which fixes the position at which the water surface is accurately 
horizontal right up to the surface of the sphere. ‘Thus, if U is the observed 
upward force, U, the upthrust for complete immersion, y the height of the centre 
of the sphere above the plane liquid surface and a the radius of the sphere, the 


necessary criterion is 

Une Shy Li/ es 

m=5-4(2) 44 A ee | al acl gee ce (1) 
‘The angle of contact @ is then given by the relation —cos0=y/a. As has been 
pointed out previously, the method possesses all the advantages of a dynamical 
method without the necessity of a knowledge of the surface tension of the liquid. 
It is therefore more reliable than a purely visual observation of the angle of contact. 


The sphere used throughout these experiments was a steel ball-bearing, 
‘0-5 inch in diameter, which was coated by immersion in molten wax to a uniform 


| 
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thickness of 0-03 mm. The paraffin wax was a high-grade B.D.H. product 


congealing at 52°c. A fresh wax surface and clean water, freshly swept with 


/ waxed barriers, were employed for each individual determination of an angle 
of contact. The position of the water surface was determined to an accuracy 


of 0-01 mm. by means of a micrometer needle, so arranged that when contact was 


_ established the grid of a triode valve was short-circuited to the cathode and an 


indication given by a milliammeter. The ring balance, employing a double 
mirror system, was used in conjunction with a telescope and a distant scale. 
A reading of 4-1 cm. corresponded to a force of 1 gm. wt., while the scale could 
be read to one-tenth of a millimetre at the velocities employed. In the absence 
of vibration, individual determinations of the angle of contact of a moving surface 
were consistent to within half a degree. The velocity quoted is always the 
velocity relative to the surface of the sphere when the criterion of equation (1) 
is satisfied. 
§3. RESULTS 

In the velocity range from 0-2 mm. per minute to 1-6 mm. per minute, very 
little change took place in the advancing angle of contact, the angles observed 
at these velocities being 108° and 109° respectively. The angles of contact at 
intermediate velocities were so closely spaced that the assignment of separate 
values was not justified. ‘The advancing angles quoted for the extreme velocities, 
however, were based on a very large number of individual runs, with freshly 
prepared surfaces, and pointed unmistakably to a small but quite definite increase 
in the advancing angle with increase of velocity. 

The effect of the time of immersion on the advancing angle was studied by 
immersing the freshly prepared sphere in water for various intervals of time 
immediately prior to the run in which the angle of contact was determined. 
‘The curve (a) in Figure 2 shows that for a constant velocity of 1:6 mm. per minute 
the advancing angle decreased with increasing time of immersion, reaching the 
value 104:5° after 45 minutes. Observations were not carried beyond this 
point on account of the risk of contamination of the surfaces. 

In dealing with the receding angle of contact, care is necessary to distinguish 
between the effects of velocity and of time of immersion. At first sight it appears 
that the receding angle increases with increase of velocity. However, a lower 
velocity necessarily implies that the waxed surface is immersed in the liquid 
for a longer period of time before the critical position defined by equation (1) 
is reached. Hence it is necessary to study the effect of time of immersion alone. 
Observations were therefore made at the highest velcoity (1-6 mm. per minute) 
after immersing the waxed surface for a suitable period of time prior to the 
commencement of the run. In the case of the receding angle, this does not of 
course involve removing the sphere from the water; and the time of immersion 
is an uninterrupted interval, ending at the instant when the criterion of equation (1) 
is satisfied. 

Figure 1 shows the intersection of the calculated curve representing equation (1) 
with the experimental curve connecting U/U, and y/a. ‘The velocity of the 
receding water surface is 1-6 mm. per minute in each case, and the effect of 
increasing time of immersion is very clearly marked. ‘The receding angle of 
contact at constant velocity is shown as a function of time by curve (8) in Figure 2. 
Clearly the receding angle decreases very considerably with increasing time of 
immersion, apparently approaching a limiting value, whose exact observation 


i 


would be difficult, however, on account of the risk of contamination of the surfaces. 
during intervals of the order of an hour or more. 
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The curve (c) in Figure 2 shows the observed receding angles corresponding» 


to different velocities as a function of the time of immersion which is unavoidably | 


associated with the measurement at that velocity. The effect of velocity alone 
on the receding angle, after particular times of immersion, may be inferred from 
the difference in ordinate of the two curves (b) and (c). Thus, when the time: 
of immersion is 38 minutes, an increase in velocity from 0-1 mm. per minute 
to 1:6mm. per minute is accompanied by a decrease in the receding angle of 
contact from 91-5° to 87:5°. At lower times of immersion the effect of a change 
of velocity is less pronounced ; and it appears that the greater part of the observed 
hysteresis is due directly to changes accompanying prolonged immersion. For 
this reason it appears justifiable to continue curve (4), as indicated by the broken. 
section, up to a receding angle not much less than the corresponding advancing 
angle of 109° for ‘‘zero”’ time of immersion. 


§4. DISCUSSION 
Comparison with the experiments of Ablett (1923) using the rotating cylinder 
method reveals a general agreement as regards the change of the advancing and 


‘receding angles with increase of velocity. Ablett was able to make observations. 


over a very much wider range of velocity, whereas the present work is confined 


of necessity to very low velocities. Ablett, however, takes no account of the 


effect of the time of immersion, though it is clear now that the effect must have 


been appreciable, at any rate at the lower end of his velocity range. Itis apparent 


also that what Ablett regards as an “‘ equilibrium” angle of contact, and determines. 
during the slow rise of the water past a stationary cylinder, is in reality an advancing 


angle. In any case, exact numerical agreement with Ablett’s work is not to be — 


expected, if only for the reason that his wax surface was prepared by turning in 
a lathe and ours by dipping in the melted wax. Moreover, Ablett does not 
appear to have used the sweeping technique in order to obtain a clean water 
surface. 4 


It is clear from the experiments described in this paper that a definite value 


can be assigned to the angle of contact of a liquid only if the exact conditions 
of observation are also specified. In quoting an angle of contact, it is necessary 
to specify not only the state of the surfaces but also their relative velocity and the 
time during which they have been in contact. There does not appear to be a 
single angle of contact corresponding to rest conditions. On the contrary, 
for relative velocities less than 0-1 mm. per minute in either direction, angles 
of contact have been observed over the range 91° to 108°. Thus, it appears. 
that the term ‘equilibrium angle of contact” is practically meaningless. Some 
theoretical interpretation of the variation of the angle of contact is attempted 
in a separate paper (Yarnold and Mason 1949). 
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* ABSTRACT. Measurements of the absorption of ultrasonic waves in water and ethyl 
} alcohol were made by the pulse method, and are more accurate than those made by other 
methods. ‘The absorption coefficient a in water was found to vary as the square of the 
| frequency v over a wide range of temperatures. "The measurements in water were made at 
) frequencies between 7:5 and 67:5 Mc/s. and temperatures between 0 and 95°c. The 
absorption decreased with increasing temperature by a factor of 8 between freezing and 
boiling points. The observed absorption was about three times that calculated from Stokes’ 
) formula and this ratio was nearly independent of temperature. In ethyl alcohol measure- 
f “ments were made at 52-4 Mc/s., and at temperatures between —50° and +60°c. The 
. absorption decreased by a factor of 3 in this range and observed values were about twice 
» those given by Stokes’ formula. Values of (a/v?) at 25° c. are 22-0 x 10-17 sec?/cm. for water 
and 50:5 x 10-1” sec?/cm. for ethyl alcohol. 

Published values of the absorption coefficients of a wide range of pure liquids and of 
some mixtures are reviewed, and a classification proposed. The experimental results 
suggest that the absorption in associated liquids is due to a different mechanism from that 
in non-associated liquids. In the latter the absorption is probably due to the same pheno- 
menon as in gases: a slow exchange of energy between different degrees of freedom. The 
relaxation frequencies must be above 200 Mc/s. for most liquids. A simple theory is given 
to explain the observed variation of the absorption with concentration in mixtures of non- 
associated liquids. The semi-quantitative agreement of theory with experiment lends 
| support to the relaxation hypothesis. 


§1. INTRODUCTION 


HE purpose of the present article is twofold. In the first place new measure- 

ments of the absorption of ultrasonic waves in water and ethyl alcohol 

are presented; use of the pulse method enabled the observations to be 
more accurate and cover a much wider temperature range than previous measure- 
ments. In the second place the results obtained by other workers in many 
different liquids are summarized and discussed. An attempt is then made to 
classify liquids by the magnitude of the coefficient of absorption and the way 
in which it varies with frequency and temperature. ‘The significance of the 
classification is discussed in relation to the molecular species of the liquids and 
it is shown that a coherent interpretation can be given to what has up to now 
seemed a diverse collection of data. A simple argument based on the relaxation 
theory is given to explain the observed variation of absorption with concentration 
in certain mixtures of liquids. 

Both the pressure and amplitude of the sound have been used in the past 
| to define the coefficient of absorption «; in what follows « is defined by the relation 
| p=p)e**, where p is the amplitude of the excess pressure at a distance « from a 
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source of plane waves of amplitude py. A quantity having perhaps more physical 
significance than « is the absorption per wavelength, ., defined by the relation | 
y.=ad. The quantity » has the dimensions of a pure number and, provided 
it is small, is also given by »=4AE/E, where E is the flux of energy per second 
through unit area normal to the direction of propagation and AF the energy per 
second per unit area absorbed in one wavelength. Measurements of the absorp- 
tion in gases acquire a greater significance if expressed in terms of yw rather than 
x, and it is considered that the same may well be true for certain classes of liquid. 

Since the development of the acoustic interferometer by Pierce (1925), 
numerous measurements of the velocity and absorption of sound in gases have 
been published; these have been summarized by Richards (1939). ‘The results 
showed that the theoretical treatments of Stokes (1845) and Kirchhoff (1868) 
were inadequate, and that the measured absorption was considerably greater 
than that predicted. 

This anomaly was treated theoretically by Herzfeld and Rice (1928), Bourgin 
(1929-1936), Kneser (1931 and 1933) and others, who showed that a large 
attenuation additional to the Stokes/Kirchhoff value was to be expected, if a time 
lag occurred in the sharing of anergy between the various degrees of freedom 
in agas. This‘‘relaxation time” 7 was attributed to the sharing of molecular 
kinetic energy between vibrational motion on the one hand, and translational 
and rotational motions on the other. An important consequence of the theory 
is that » should increase to a maximum value at a frequency »,,, related to the 
relaxation time by t=1/27v,,.. The theory also predicts a change in velocity in 
the same range of frequencies. ‘These two effects have been well substantiated 
by experiment for polyatomic gases (Richards 1939). 

The first measurements of velocity in liquids were those of Freyer, Hubbard 
and Andrews (1929), and the first reliable work on absorption was that of Biquard 
(1935). Other data on absorption published in the thirties (cf. Richards 1939) 
are often not reliable, as indicated by the disagreement between different authors. 
No conclusive evidence of a change in velocity with frequency was found, but 
all workers were agreed that the absorption was often much greater than that 
to be expected from viscosity according to Stokes’ formula; in some liquids the 
observed value of « was 100 times that calculated. (In most liquids the attenuation 
calculated from thermal conductivity (Kirchhoff 1868) is only a small fraction 
of that due to viscosity; mercury and liquid helium are exceptions to this). 

Many attempts have been made to explain this anomaly for liquids in the 
same way as for gases, but the absence of the two effects of a relaxation phenomenon 
mentioned above appeared to be a serious obstacle to the theory. 

Kneser (1938) showed however that published observations could be 
explained by postulating a relaxation phenomenon of a general nature, with a 
relaxation frequency v,, so far above the experimental frequencies that the 
expected dispersion effects would not have been observed. 

The problem can now be re-examined in the light of the more accurate experi- 
ments which have become available (Willard 1941, Pellam and Galt 1946, Pellam 
and Squire 1947, Pinkerton 1947). The author’s measurements in water and 
ethyl alcohol described in §2 yield results differing from those obtained by other 
workers for non-associated organic liquids like carbon disulphide and benzene, 
in that the absorption is much nedter to the Stokes value and the temperature 
coefficient is larger and negative. 
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§2. EXPERIMENTAL METHOD AND RESULTS 
(i) Method 


‘The experimental arrangement and technique will be described in a separate 
/ communication, and it will be sufficient here to indicate the general principle 
| ofthe method. A short account of the measurements in water has already appeared 
| (Pinkerton 1947). 

. The pulse method, first described by Pellam and Galt (1946) was developed 
) independently on similar lines by the author from the ultrasonic techniques 
4 used for radar during the war. In this method short bursts of ultrasonic oscilla- 
) tions are excited in a quartz crystal. Pulses of ultrasonic waves travel out into 
| the liquid from one face of the crystal, and are returned by a movable reflector. 
1 The reflected ‘“‘echo”’ is reconverted into electrical oscillations by the same 
| crystal and displayed, after amplification, as a deflection of the trace of a cathode 
» ray oscillograph. Measurement of the amplitude of the echo as a function of 
' range yields a value for «. Due attention is paid to the diffraction effects, which 
_ cause a spreading of the beam from the source, and special care is taken to control 
the temperature of the liquid. The pulse method has the advantage, not possessed 
| by all other methods, that the power transmitted is so small that errors due to 
) local heating and cavitation are negligible. 
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Measurements were made over a wide range of temperatures, in water from 
0 to +95°c., and in alcohol. from —50° to +70°c. The apparatus was also 
designed so that measurements could be made on several frequencies with the 
same crystal by the use of harmonics. The frequencies used were approximately 
7-5, 22:5, 37-5, 52-5 and 67:5 Mc/s.; a second crystal was used for measurements 
at 15 Mc/s. Pulses of duration varying from 2 to 20 microseconds were used. 


(ii) Results 


| (a) Water. Measurements were made first on 15-0 Mc/s. ‘To obtain a 
_ path corresponding to a total absorption of 40 db., a tank 150 cm. long was thought 
desirable. Unfortunately the losses of heat from so large a tank made it imprac- 
ticable to observe above 68° c. with the thermostat available. ‘The corrected 
results at 15 Mc/s. are shown in Figure 1. 
In order to find the variation of « with frequency v measurements were later 
made at the other five frequencies quoted, at temperatures between 0° c. and 
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about 60° c. and curves similar to Figure 1 were drawn for each frequency. 
These curves were used to find the variation of « with v at the temperature up 


to 60° c. quoted in Table 1; Figure 2 shows values of « as a function of v at three © 
different temperatures. The frequency was not always an exact harmonic of | 


7:5 Mc/s. but was measured separately during each experiment with an absorption 
wavemeter incorporating a crystal check circuit. 

It can be seen that « varies rapidly with v and a plot of «/v? against v, shown 
in Figure 3, indicates that «/v? is approximately constant at constant temperature. 
The accuracy to which «/v? remains constant can be assessed by using the depar- 
tures of [«/v?], from the mean values to calculate the probable error. The 
results of these calculations are shown in Table 1 and justify the statement that 
a/v? is constant to the order of accuracy of our experiments. 

For the measurements of « between 60° c. and 95° c. a smaller trough was 
used to reduce thermal losses, and the shorter path length available restricted 
measurements to the single frequency of 52:5 Mc/s. 

The most striking feature of these results is seen by comparing the observed 
value of «/v? with that predicted by the theory of Stokes (1845), who showed 
that * «/v?=877n/3pc? where 7»=coefficient of shear viscosity, p=density, 


Table 1. Absorption of Sound in Water 


Measured values compared with those predicted from Stokes’ formula 


(1) (2) (3) (4) (5) 
T (a/¥°)ot (a/>*)eat a) ee 
(COP aGec7 cma 0m) (Gees /em= a liias4) Seals Ay 
Qobs (%) 
0) 56-9 17-09 30338! 0-6 
5 44-1 13-40 3°29 0-6 
10 35°8 11-34 sya) 0-86 
15 29°8 9-59 )0H | 0-9 
20 DSS 8-16 3-10 1-35 
30 19-1 6-15 3-10 1:7 
40 ; 14-61 4-85 3-01 0:5 
50 11-99 3:98 3-01 0-4 
60 10°15 3-36 3-00 0-8 
70 8°71 2:91 2-99 = 
80 7:89 ’ 2:58 3-06 — 
90 7°24 ASSIA 3-14 — 
100 6:87 214 3:26 — 


Values of a/v? up to 60° c. are means of measurements at six frequencies in the range 
7:5 to 67:5 Mc/s.; values from 60 to 100 are from measurements at 52°5 Mc/s. only. An 
extrapolation to 100° c. has been made from 95° c 


c=velocity of sound. Column (3) of Table 1 gives the values of «/v? calculated — 


from this formula for various temperatures and column (4) the ratio of observed 
to calculated absorption. ‘This ratio is strikingly constant over the whole range 
of temperatures. Values of 7 and p used to find («/v?).,;. were taken from the 
International Critical Tables, and the velocity from Randall (1932) whose 
measurements are probably the best for water. 

* This is actually a simplified form of Stokes’ result obtained by assuming, on grounds elaborated 
by him, that any viscosity effects involved in a pure compression could be ignored. The present 


results may be interpreted as an indication that the second or compressional viscosity is not zero in 
water or alcohol. 
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(b) Ethyl Alcohol. Since Pellam and Galt (1946) have stated that the 
absorption in ethyl alcohol is less than is calculated from the Stokes formula, 
| measurements were also made in this liquid. Two samples of alcohol were 
2 used; the first contained 6°, water and the second not more than 0:5%, and 
| they were investigated at a single frequency of 52-4 Mc/s. The results of the 
| measurements on the dry sample are given in Table 2 in which the values of 
| (a/v?) op; have been taken from the best smooth curve through the experimental 
points. The absorption in wet alcohol was the same as in the dry sample at room 
) temperature but rose to a value 30% higher at —30° c.; at +60° c. on the other 
| hand it dropped 5% below that of the dry sample. 

The results for the dry sample are compared with Stokes theory in column (7) 
{ of Table 2. The velocities of sound for the lower temperatures were obtained 
_ by linear extrapolation of the results of Freyer, Hubbard and Andrews (1929), 
§ and the values of viscosity interpolated from the results of ‘Tonomura (1933). 
) Above 0° c. the data for 7 and p were taken from the International Critical Tables. 


Table 2. Absorption in Ethyl Alcohol 


) (2) (3) (4) (5) (6) (7) 
r (a/V*)ovs Uf] Velocity (a/V*)eale obs 
(rc) (sec?/em=< 101") (poise) p (cm/sec. X 10°) (sec?/em. x 107") quar 
—50 137-0 0:0630 0-8497 1-418* 68-4 2-00 
—40 113-2 0:0475  —0-8405 ogres 56-3 2-01 
— 30 93-2 0:0367 0:8319 1:346* 47-6 1:94 
—20 78-2 0-0287 08233 fhegititis 40-7 1-92 
—10 67-6 0-0226 0:8147 (ears 32 1-92 
0) 60-4 0-:01790 0-8062 1-242 30-5 1:98 
10 Sy) 0:01451 0-7978 1-204 27-4 2:04 
20 52-0 001194 0-7893 1-168 24-96 2-08 
30 48-7 0-00992 0-7810 1-134 22-90 2-18 
40 46:5 0:00831 0:7727 1-101 Pe) 2:19 
50 44-9 0:00697 0-7633 1-067 19-79 227 
60 43-7 0:00596 = 0-7561 1-034* 18-76 D233 
70 42-4 0:00504 0-7470 1-000* 7-8 2:39 


* Extrapolated values. 


‘Table 2 shows that the observed absorption is approximately twice the calculated 
at all temperatures. This result does not agree with that of Pellam and Galt 
who found absorption less than the calculated value; it agrees, however, with 
the observations of Willis (1947) within his experimental error. Since the 
variation of « over the range of temperatures used is 8:1 in water and only about 
3:1 for alcohol the constancy of «). /%eaj. may be considered more remarkable 
in water. 


cale 


§3. REVIEW OF PUBLISHED MEASUREMENTS OF ABSORPTION 
(i) Experimental Results for Pure Liquids 


To establish that the absorptions in different liquids are in accord with certain 
general rules, we now survey the measurements published by other workers, 
whose values for a variety of liquids are collected in Table 4. Values of «/v? 
and not « are given and form a fair basis for comparison, since it has been shown 
for most of the liquids in the table that «/v? is independent of v (Bar 1937, Grobe 
1938, Willard 1941, Rapuano 1947 and others). A decrease in «/v? with 
ncreasing v seems indicated by the results of Rapuano (1947) for carbon 
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disulphide; at 75 and 105 Mc/s. «/v? was about 1/5 of its value near 10 Mc/s. | 
(Biquard 1935, Willard 1941). -«/v® is far from constant for acetic acid, methyl — 


and ethyl acetate and solutions of metallic acetates (Bazulin 1936, 1939, Claeys, 


Errera and Sack 1937, Pinkerton 1948); a full discussion of these latter cases is f 


not relevant here and will therefore be reserved for a future communication. 


With the above exceptions it may be taken as established that « is accurately — 


proportional to v? in most liquids up to 100 or 200 Mc/s. At higher frequencies 
evidence is lacking except for mercury for which «/v? is constant up to 1,000 Mc/s. 
(Ringo, Fitzgerald and Hurdle 1947). 

The values of Table 4 are taken from Willard (1941) and Pellam and Galt 
(1946), who measured in all about 50 organic liquids; their results agree within. 
the experimental accuracy claimed, and are confirmed for the commoner liquids. 
by Biquard and others. Results are also available for the two monatomic 
liquids, helium (Pellam and Squire 1947) and mercury (Bar 1937, Rieckmann 
1939, and Ringo, Fitzgerald and Hurdle 1947), in which the absorption was 
found to agree with the Stokes/Kirchhoff value (helium shows an anomaly near and 
belowtheA-point). Four highly viscous liquids were investigated by Hunter (1941). 

The dependence of « on temperature is highly significant for any classification 
of liquids according to their absorption. This dependence was investigated 
for many liquids by Pellam and Galt (1946) and also by other workers for a 


number of individual liquids, viz.: water and alcohol by the author (v. supra), 
water (Fox and Rock 1946), benzene (Quinn 1946, Bazulin 1937), carbon 
tetrachloride (Bazulin 1937), four very viscous liquids (Hunter 1941) and liquid — 


helium (Pellam and Squire 1947). 


(ii) Classification of Liquids by their Absorption 
A previous attempt has been made by Kittel (1946) to classify liquids according 
to their absorption of ultrasonic waves. A more extensive classification is here 
proposed based on two criteria: (a) the ratio of observed absorption to the 
Stokes value, denoted «}./%:a;.; (6) the temperature coefficient of absorption. 
It is found that, with the exceptions mentioned in §3 (i), all pure liquids may be 
assigned to one of the four classes of Table 3. 


Table 3 
(1) (2) - 3) (4) (5) 
Z Range of ‘Temperature Type of 
cee (Gobs/@cale) coefficient of a liquid Examples 


Anomalous 1,500 to 3 Positive. (aops/@eate) Non-associated Carbon disulphide; 
Al varies with 7. polyatomic. benzene ; most 
non-associated of- 
ganic liquids. 


Anomalous Stone Strongly negative. Associated Water; alcohols. 
All (Qops/Geale) nearly polyatomic. 
independent of T. 
Normal i Positive for helium. Monatomic. Helium; mercury. 
NI 
Normal i Negative. Associated Glycerine ; highly 
N {I polyatomic. viscous liquids. 


‘Table 4 shows a number of liquids listed according to the classification of Table 3 ; 
in each class the liquids are quoted in descending order of a/v. 
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The two anomalous classes of liquids accord well with the more usual dis-_ 
tinction between associated and non-associated liquids. Moreover the liquids 
in class Al with the greatest absorption are those having zero dipole moment 
and small roughly symmetrical molecules. The introduction of chlorine into 
the benzene ring to form chlorbenzene destroys the symmetry and lowers the | 
absorption whilst increasing the dipole moment. ‘This tendency is more marked 
in nitrobenzene whose absorption is still less. In class A II the nearly regular 
progression of values of «,),, from water to amyl alcohol is most striking, although 
butyl alcohol does not seem to take its proper place in the series, probably because 
the absorption quoted for this alcohol relates to 2° c. and not to some temperature 
near 25° c. Even more noteworthy is the fall in the ratio «,),./%cai, aS one ascends 
the homologous series, water taking the place of an “alcohol” with no carbon 
atom. It has been mentioned that in water and ethyl alcohol «,),./%.9). is nearly 
independent of temperature; sufficient data are not available to test this rule 
for the other alcohols. 


(iii) Experimental Behaviour of Mixtures 


The absorption of ultrasonic waves in mixtures of two liquids in varying: 
proportions has been investigated in the following cases: mixtures of benzene 
with toluene (Bazulin 1936, Grobe 1938), carbon disulphide with benzene and 
with carbon tetrachloride (Claeys, Errera and Sack 1937), and benzene with 
acetone (Willard 1941). All these liquids are of class AJ and are non-associated ; 
acetone and toluene have dipole moments. In each case one of the pure liquids 
had an absorption about ten times larger than the other. 

In all the above cases there was a very rapid decrease in « as the proportion of 
the less absorbing liquid was increased (cf. Figure 5); further increases in the 
proportion of the latter had progressively less effect. ‘There was no suggestion 
of a linear variation with concentration. 

Contrasting results were obtained by Bazulin and Merson (1939), Willard 
(1941) and Willis (1947), who measured the absorption in mixtures of acetone - 
and water, and by Mikhailov (1940) Mikhailov and Gurevich (1946), and Willis 
(1947) who studied mixtures of water and ethyl alcohol. Mikhailov and Gurevich 
also studied mixtures of methyl alcohol and water. All these observers found a 
maximum in the absorption at some intermediate concentration, the maximum 
absorption being about twice that of the more absorbing liquid when pure. 


<4. THEORETICAL DISCUSSION 
(i) Possibility of Scattering 
It has occasionally been suggested (e.g. Lucas and Biquard 1937) that the 
anomalous part of the absorption is due to some scattering process. It is a simple 
“matter to test for scattering using the pulse technique, and special tests were 
made in water and alcohol free from dust and air bubbles. ‘These tests failed 
to reveal any scattered sound, although if scattering were to account for any 
appreciable part of the absorption it should readily have been perceptible under 


the conditions of the experiment. We conclude that the anomalous values 
of the absorption coefficient are not due to scattering. 


(ii) The Relaxation Hypothesis 


The relaxation phenomena found in gases and caused in that case by the slow 
equipartition of energy between different degrees of freedom, will also occur 
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| in any other case where the volume changes of the medium do not follow instant- 
| aneously the changes of pressure in the sound wave. Kneser (1938) showed that 
| ifa single relaxation process gives rise to absorption in a medium, then the expected 
| variation of « with v would be of the form 


re at yee he ee (1) 


where A is very nearly constant, and v,,=1/2m7, is the postulated relaxation 
| frequency. From (1) «/v? would be expected to be constant so long as v<v,, 
and from the experimental values of «/v? Kneser made estimates of the lowest 
, possible relaxation frequency in various liquids. With the data now available 
) his estimates of v,, would have to be increased in many cases since we know «/v? 
to be constant up to 100 or 200 Mc/s. in several liquids (Rapuano 1947). 

However we may set an approximate upper limit to v,, as follows; we consider 

the absorption per wavelength, u, which from (1) is given by 


f= aA=(a/v®)®A= Acr|[1 + (o/%)°] = BO) Flr 2] eee (2) 


where B is a new constant = Acy,,, and c is the velocity. Equation (2) shows that, 
for v<v,,, wu increases linearly with v as on the Stokes theory. When v=»,, 
however mw attains a maximum value, p,,=Acv,,/2. Since A and ¢ are known 
from experiment, an éstimate of z,, gives a value for v,,. It is possible to make a 
reasonable guess at y,, for a liquid if the absorption is due to slowness of exchanges 
of energy between degrees of freedom, from the value obtained for ,, in gases. 
If this mechanism be accepted p,, is completely determined by the specific 
heats C,, and C, and the fraction of C, corresponding to internal degrees of 
freedom (Richards 1939). Since these specific heats probably do not change 
by large factors on passing from gas to liquid (Herzfeld 1941), we put 
(Hm) gas (Mmiquae Values of (2u,,) quoted by Fricke (1941) vary between 0-2 
and 0-3 for several gases, and in other cases can be calculated from specific 
heats of the individual degrees of freedom obtained from spectroscopic data .* 
'To find »,, we take 2u,,~0-25, and therefore from (2), Acv,,==0-25 or 


POAC eR ey a (3) 


Here A is the constant value of «/v? found in experiments at low frequencies ; 


tentative values of v,, calculated from (3) are given in Table 5. 


m 


Table 5. Approximate Values of v,, for 2u,, 0-25 


Liquid v,, (Me/s.) Liquid Vy, (Me/s.) 
Carbon disulphide 30 Nitrobenzene 2,200 
Benzene 240 Ethyl] alcohol 4,000 
Carbon tetrachloride 500 Methy] alcohol 6,000 
Chlorbenzene 1,500 Water 7,000 
~ "Toluene 2,000 


‘The value v,,=240 Mc/s. shown for benzene is sufficiently small to deserve 
comment; measurements were made up to 165 Mc/s. by Rapuano (1947), and 
in the diagram in which his values are given there seems to be a small but systematic 
decrease in «/v? up to 165 Mc/s., which does not conflict with our estimate for »,,. 
It is more surprising that the relaxation frequency of 30 Mc/s. predicted 

for carbon disulphide has not been found experimentally; the drop in «/v? 


* Unpublished calculations for liquids by Mr. E. Bauer from spectroscopic data have also shown 
that this is a reasonable estimate for pm. 


= 


138 F. M. M. Pinkerton 


between 10 and 75 Mc/s. mentioned in §3 (i) is partial confirmation, and a more 
detailed investigation is planned by the author. 

The experimental results have therefore been shown to be consistent with 
the occurrence of a relaxation process at a frequency outside the present range 
of measurement. The following facts indicate that the relaxation mechanisms. 
are different for the liquids in classes AI and ALI. 

(i) The temperature coefficient of absorption is positive for class AI and 
negative for AIL. In class AIL «,),./a¢ai¢ 18 nearly independent of 7. 

(ii) The absorption varies very differently in mixtures of liquids in the two. 
classes. 

(iii) Calculations by Herzfeld (1941) showed that the specific heat of the internal 
vibrations in water is too small to account for the anomaly and that the absorption — 


Pf ale em amie = 


from this cause should be zero at 4° c. when C,,=C,, which is not true (Fox | 


and Rock 1946). , 

It is suggested that relaxation of the internal degrees of freedom postulated 
for gases will account for the anomalous absorption in liquids of class AI. The 
behaviour of liquids in class AIL is probably determined by changes in local 
molecular order (Hall 1947, 1948). In certain other liquids like acetic acid 
the absorption may be due to the perturbation of a chemical equilibrrum 
(Pinkerton 1948). ’ 

A complete theory of the absorption of sound in liquids in class AI is required 
to predict v,, as well as y,,, but there seems little prospect at present of being 
able to predict v,, from independent data. However the theory outlined above 
implies that the wide variations in «/v? between liquids must be largely due to. 
variation of v,,, since we have assumed j1,, to be roughly the same in all liquids.. 
From equation (2) it is easily seen that 


Lim (a/y)=2n ev.) ON) eee (4) 


v—> 0 


This relation shows that highly absorbing liquids must have small values of vy, — 


which implies that they do not readily share their vibrational energy. But this. 
is what we should expect for liquids having small spherically symmetrical mole- 
cules with no dipole moment ; such liquids should therefore be strongly absorbing. 
Reference to class AI of ‘Table 4 shows this to be the case. The symmetrical 
non-polar molecules of carbon disulphide, benzene and carbon tetrachloride 
show the greatest absorption, while nitrobenzene and propyl iodide show small 
absorption. 

The same general rule was followed in many other liquids omitted from 
‘Table 4 for reasons of space. 

The frequency of molecular collisions in a liquid is about 10" per second 
at room temperature (Glasstone, Laidler and Eyring 1941), so that taking the 
relaxation frequencies to be of the order of those in Table 5, it follows that from 
500 to 30,000 collisions are required before the internal and external energies. 
come to equilibrium after an abrupt change in the external conditions. It is 
noteworthy that in those gases, for which the lifetimes of the excited states are 
longest, v,, and therefore the absorption are most sensitive to traces of impurities. 
(Fricke and Knudsen 1941), although y,, is not greatly affected. It has frequently 
been suggested that this is due to the impurity acting as a “‘ catalyst” in transfers 
of internal energy. Such a process may very well also occur in mixtures of two 
liquids, and is the basis for the argument used below. 


——. 
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(11) Explanation of the Behaviour of Mixtures 


By making certain simple assumptions the relaxation mechanism proposed 
can be used to calculate approximately the variation of absorption in mixtures 
| of two liquids of class AI. We introduce a relaxation time 7 =1/27v,, where: 
| Vy, 18 the relaxation frequency of equations (1)-(4). In a mixture of liquids A 
| and B, in which the equilibrium between the energies associated with internal 
and external degrees of freedom is set up by collisions, it is possible to define 
four relaxation times as follows: 7, , and rx» are the relaxation times for collisions. 
between similar molecules of types A and B respectively, +, is that for collisions. 
of excited type A with de-excited B, and 7,, that for collisions of excited B with 
de-excited A. Since the absorption in A is much greater than in B we now 
BssUIle-7,, -Typ and Write 7,,=Tp,=Tpy. Lhe reason for this is that if 
molecules of type B possess a stronger mutual interaction, the interaction between 
A and B molecules may well be of the same order as between two molecules of 
B, and may be expected to be much greater than the mutual interactions between 
molecules of type A. Denote the fraction of molecules of type B in the mixture: 
by x, and neglect differences in the molecular diameters, which is justifiable 
to a first approximation for the liquids on which data are available. Then in 
unit time a fraction (1 —x) of A will collide with molecules of the same kind and 
have relaxation time 7,,, and a fraction x of A will have relaxation time T,p,. 
because of collisions of A with B. Ina short time Az the number of molecules A 
relaxing will be proportional to: 


At At At 
—(l-—x)+—x=—, say,  —— ...... (6) 
TAA TBB TA or 

where 7, ,. is the effective relaxation time for molecules type A. Now from the 


assumption that ju,, is the same in both pure liquids, it follows from equation (4): 
that, at low frequencies 


t,/%_= (CBY mg) /(Ca Vm) = (Cp744)/(CaTBp)- 


It will be convenient to assume that c, =cp, which is true for the case of benzene 
and toluene. We can then put 

Lyn — Ti Tap = 1/y, Say. OS oOo (7) 
In the mixture the absorption will be additive for the two molecular species and. 
will be given from equation (4) by: 


on(x) ue ie [75 Al —x) +7Tppx]. ireereo25 (8) 


py 
Substituting for 7, _. from (6), and using the conditions that 
o(x)=o,., when «=0; axj=an, when v=], 
we arrive at the result that: 
| 
a EX, Vong SAV. oso tr 2 
o(s)=y] ESA +0 | 24 = Ply) ty say.) 


This expression gives the proportional variation in absorption in the mixture, 
in terms of the ratio 1/y of the absorptions in the pure liquids, and the molecular 
proportion x. Both quantities are known from experiment, and no arbitrary 
parameter need be introduced. 
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In Figure 4 the relation between P and x is shown for the two arbitrary values 

of y, 0:3 and 0:1. A detailed comparison with experiment is made in Figure 5 

in which the points indicate values found by Grobe (1938) for mixtures of benzene 

with toluene. The full curve corresponds to y=0-1, which is the ratio of the 
absorptions in the pure liquids. 


Fraction of Meximum Absorption 


06 08 10 


02.04 j (020-4 wae Osean 
Molecular Proportion 2 ‘ Molecular Proportion of Toluene 
Figure 4. Graph oe function Figure 5. Observed and theoreti- 
x cal absorption in mixtures 
= ee ee, 
PX, y)=y (1—x)y+x | x| of benzene and toluene. 
far 02a f and y=071 oe curve primey and 
meee the points due to Grobe 
(1938). 


The agreement between observed and calculated values is satisfactory in 
view of the drastic simplifying assumptions made. Equally good agreement 
is found for the other pairs of unassociated liquids which have been investigated. 

This explanation does not fit the case of mixtures of associated liquids at all; 
the process which is responsible for their behaviour is not clear at present. 
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DEAR S LO THE EDITOR 


A Note on Magnetic Viscosity in Alnico 


Magnetic viscosity occurring in specimens of a high-coercivity alloy, alnico, has already 
been briefly reported (Street and Woolley 1948). It would appear from further investi- 
gations carried out in this laboratory that the experimental results may be satisfactorily 
accounted for by assuming that, after a sudden change in an external magnetic field 
applied to the specimen, the movement of the domain magnetization vectors is restricted by 
energy barriers. Thus the domains may be characterized by a distribution of activation 
energies, and by using an analysis similar to that of C. L. Smith (1948) for transient creep 
in metals we find that the total change in the intensity of magnetization, AJ, at a time 
t sec. after the change in field has been completed is given by the equation 


AT=(constant)T log (C2), 


where T is the absolute temperature of the specimen and C 1s a constant. 
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Figure 1. Results for alnico for Figure 2. Ewing’s results (1889) for 
H = coercive field. soft iron. 


It would be expected from this equation that AJ would be a linear function of log f, 
and the (AJ, log t) curve obtained with an alnico specimen at its coercive force point 
(Figure 1) indicates that the experimental results are in satisfactory agreement with this 
hypothesis over a large range of t. The results of Ewing (1889) for soft iron specimens are 
represented in the form of a (AJ, log t) curve (Figure 2), and here also a linear variation exists. 
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Again, the equation predicts that the slope of the (AJ, log t) curve should be directly 
proportional to the absolute temperature T of the specimen. Preliminary experiments with” 
an alnico rod seem to support this conclusion. 

It is hoped to publish a more detailed account of this work at a later date. 


The University, Nottingham. R. STREET. 
25th November 1948. J. C. Woo.tey. 


Ewine, J. A., 1889, Proc. Roy. Soc. A, 46, 269. 
Situ, C. L., 1948, Proc. Phys. Soc., 61, 201. 
SrreeT, R., and Woo.tey, J. C., 1948, Proc. Phys. Soc., 61, 391. 


REVIEWS OF BOOKS 


Shadow and Diffusion in Illuminating Engineering, by K. NoRDEN. Pp. x+165. 
(London: Pitman, 1948.) 25s. 


The quality of lighting denoted by hardness or softness is a very familiar one. On it 
depend the sharpness or vagueness of shadows, and the contrasts and prominence of shading 
which help to reveal the shapes of objects seen. The very great importance of the correct 
disposition of shadows, and of correct contrasts, is seen in the work of a skilled photographer, 
‘or in considering the conditions necessary for the best visibility in various industrial tasks. 
The recent work of H. C. Weston is notable as having established the relation between 
contrast and similar factors on the one hand and ease and speed of performance of visual 
tasks on the other. It has done much to increase the interest of the lighting engineer in 
the whole subject of the directional properties of lighting systems in relation to the appear- 
ance of the objects illuminated. It is not a simple matter, however, to translate the degree of 
‘contrast and other factors governing the visibility of an object or task into corresponding 
factors in the design of a lighting installation as a whole, suitable for inclusion in the 
engineer’s specification. Moreover, the lighting designer has in any case many other factors 
to think about, such as the level, diversity and colour of the illumination, and glare. He may 
well claim, therefore, that to try to measure or specify in detail the vectorial distribution of 
light on the working plane, in any terms but those of a general commonsense empiricism, 
would be really too much of a good thing. e 

Much of Dr. Norden’s book is an open challenge to this attitude. Dr. Norden proposes 
amore formal approach, through the study and measurement, according to certain standard 
methods, of the shadow-casting properties of the lighting. 'The proposed methods have 
the merit of directness and of affording useful information of the kind required without 
unduly complicated analysis. On the other hand it remains to be proved that the use of the 
methods would lead to better results than those now achieved by a simpler approach, and 
also that the lighting engineer would find the methods sufficiently manageable and profitable, 
to outweigh the extra complexity of an already difficult task. One has considerable doubt, 
therefore, whether the challenge will be accepted. 

The rest of the book is concerned, broadly speaking, with two other aspects of the subject. 
"The first is the nature and formation of shadows. ‘This is treated in the main geometrically, 
and many of the examples and calculations given are either self-evident or not original. 
‘The second is the consideration in more general terms of the revealing power of light and 
‘shade, and their relation to the visibility of fine detail, and so on. It is in the treatment of 
this aspect that one feels Dr. Norden really arouses one’s interest and enthusiasm, and 
‘succeeds in his advocacy of the importance of his subject. 

Dr. Norden has been known to lighting engineers for many years as a courageous and 
dogged protagonist of the study of shadows. ‘The fundamental rightness of his course could 
not now be denied; but on reading this latest publication of his one cannot help wishing that 
it had been devoted only and entirely to a study of the subject for its own sake, and not with 
a view to providing new tools for the installations designer. In treating this broader aspect 
of his subject Dr. Norden combines the analytical approach of the scientist with the sensi- 
bility of the artist. A book devoted simply to examples of visual tasks drawn from industry 
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and common experience, well illustrated with plates like the excellent though not numerous 
‘ones in Dr. Norden’s book, and with commentaries by an author having both scientific and 


artistic ability, would be most welcome. 


Finally, it is not out of place to say that one looks forward to the day when the attractive- 


ness of the binding of a book such as the one here reviewed can be matched by the quality 
of the paper available for the text. J. S. PRESTON. 


Theoretical Structural Metallurgy, by A. H. CorTrrety. Pp. vii +256. 
(E. Arnold & Co., 1948.) 21s. 


‘The rapid progress of the physics and mechanics of solids in recent years has placed the 
majority of metallurgists in a difficult position : with the slender furdamental knowledge 


i) dispensed in most Univeisity metallurgy degree courses (Dr. Cottrell’s University, Birming- 


ham, is a praiseworthy exception) they can hardly follow recent scientific developments 
in their own field, let alone. take an active part in them. Most of the pioneers in the 
science of metals have been physicists, chemists, or engineers by training; the ectivity of 


_ the professional metallurgist is usually confined to the application of available scientific 


results to particular metals and alloys, unless he is capable of acquiring the necessary 
fundamental knowledge by home study. In this difficult situation the appearance of 
Hume-Rothery’s brilliant non-mathematical Atomic Theory for Students of Metallurgy in 
1946 was an event of great importance, and Dr. Cottrell’s new book represents a further 
welcome addition to the means by which the metallurgist can get acquainted with the modern 
aspects of the science of metals. 

Dr. Cottrell’s book falls into two parts of rather different character. The first six 
chapters give a good elementary account of the electron theory of metals, with a brief 
sketch of the Bohr theory and of the basic conceptions of wave mechanics; the remaining 
eight chapters deal with statistical thermodynamics, heterogeneous equilibria, thermo- 
dynamical properties of metals and alloy systems, order—disorder transformations, diffusion, 
and the kinetics of phase changes. ‘These eight chapters treat their subject on a University 
science course standard; in fact, the amount of mathematics used is sometimes above the 
minimum necessary for the purpose. 

Physicists may question whether a useful purpose is served by a mainly descriptive 
treatment of a mathematical subject such as is the electron theory of metals. In the 
reviewer's opinion the answer is definitely in the affirmative. A recent text-book entitled 
Modern Metallurgy for Engineers (New York: Pitman, 1941) gives a list of references to 
literature on metals, with brief commentaries on each item, and here Mott and Jones’ book 
receives the following remark: “‘ For enthusiasts on quantum mechanics. A few engineers 
may be interested in the mental calisthenics’. Since the author of this revealing utterance 
is an influential metallurgist, and since the future of metallurgical education depends very 
much on influential metallurgists, little progress can be expected before these are convinced, 
by means of clear elementary reviews, that the physics of metals is no mere mental 
calisthenics unworthy of a he-metallurgist. 

No doubt many metallurgists will take advantage of Dr. Cottrell’s excellent book, and a 
new edition may have to be printed in due course. ‘This would give an opportunity to 
increase the accuracy of expression in many passages which, although usually correct from 
a legalistic point of view, appear likely to be misunderstood by a beginner. One of tke 
first examples of this is found on page 5, where the uncertainty principle might seem to have 
preceded quantum mechanics, in which “it has become necessary, however, to specify 
four SIE numbers instead of the single one used in the original feemaulenes of the 
Bohr theory ”’ 

On page 34 the law of rational intercepts is thought to govern the positions, not of the 
natural crystal faces, but of lattice planes (for which the statement of the law would be 
partly trivial, partly—in its emphasis on small integers—wrong). 

On page 43 we read: ‘‘ There is some evidence, particularly from the study of crystals by 
x-ray diffraction, suggesting that a real crystal is not of perfectly uniform orientation, but 
is composed of a mosaic of small crystals, or crystallites.” This seemed a possibility until 
1934, when Renninger, using the double-crystal spectrometer, found that a synthetic NaCl 
crystal was ideally perfect in the sense of the Darwin-Ewald theory; recently Guinier, 
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by means of another x-ray method, came to a substantially identical conclusion for 
recrystallized aluminium. Thus, a ‘“‘mosaic”’ imperfection is characteristic of a bad 
crystal, not of a real crystal; the early x-ray workers found most crystals imperfect because 
they did not know how to obtain perfect ones. : 

It is a pity that the last chapter, ‘‘ Phase precipitation -by nucleation and growth”’, 
does not include strain transformations (martensitic transformations), which are of the very 
greatest fundamental and practical importance in metallurgy. In this chapter Widman- 
stiitten structures are attributed to oriented nucleation and growth governed by the condition 
of minimum interface energy. That the typical Widmanstatten pattern does not arise 
by such a process can be recognized from the very frequent occurrence of lamellae penetrating 
one another; the absence of isothermal growth, the reversibility of some strain transforma- 
tions (e.g. the transformation of metastable B-brass by cooling), and many other phenomena 
show sufficiently that there can be no question of an explanation on the basis of oriented 
growth. Readers of this chapter will do well to follow it up by consulting Chapter XXII 
in Barrett’s Structure of Metals, where a balanced view of the subject is given. 


These remarks, however, should not obscure the fact that Dr. Cottrell’s book is a most | 


useful contribution to metallurgical literature; it deserves to be present in the library 
of every metallurgist, and it would certainly have deserved a more attractive paper and 
binding. E. OROWAN. 


Paramagnetic Relaxation, by C. J. Gorter. Pp. vii+127. (Amsterdam: 
Elsevier Publishing Co., 1947.) 12s. 


This monograph, by the Professor of Experimental Physics in the University of Leyden, 
was written during the last winter of the German occupation of the Netherlands, whea not 
only the absence of electricity, but even the risk of going into the streets, failed 


: 


to put a stop to the academic work of the Physics Department of the University of — 
Amsterdam. This little-known branch of magnetism, which deals with the magnetic — 


analogue of dispersion and absorption in polar substances, is treated comprehensively from 
both the experimental and theoretical standpoints, and includes much unpublished work 
carried out in Holland during the war. 
After an introductory chapter surveying the background of the work, the experimental 
methods, calling for magnetic susceptibility measurements at high frequencies, are reviewed. 
The cogent features of the different methods and their results are discussed critically in 


Jot 


the light of the new viewpoints, emphasis being placed on the various problems, both 


solved and unsolved, connected with adiabatic demagnetization and the theory of the 
crystalline state. ‘The work constitutes a remarkable achievement under the most difficult 
conditions, and can be recommended to those seriously interested in this branch of para- 
magnetism. W.S. 


Vibration and Sound, by P. M. Morse. Second Edition. Pp. xx+468. 
(New York and London: McGraw-Hill, 1948.) 33s. 


The first edition of this book earned its place among the standard theoretical text-books 
on sound. In this, the second, edition the author has taken advantage of war-time develop- 
ments in applied acoustics—to which he himself contributed not a littl—to add material, 
particularly in the subjects of radiation, diffraction and transient phenomena. ‘The latter 
especially is of great value in the study of musical instruments and auditorium acoustics. 
One is pleased to observe that sound is not treated as a watertight compartment in the vessel 
of physics, treatments being taken bodily and adapted from the theory of atomic radiation 
and of waveguides, for example. 

Not the least valuable part of the text are the frequent thought-provoking problems 
to put the student on his mettle—he will not find the going easy. Operational. calculus is 
employed in certain sections. However, with some judicious skipping, indicated by the 
author, the tyro is brought up gently from fundamental principles. Having thus scanned 
the field, it will be surprising if he does not start afresh, with his appetite whetted for a more 
detailed survey. ‘ EvGars 
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Techniques Générales du Laboratoire de Physique, Vol. 1, edited by J. SuRuGuE. 
Pp. 433. (Centre de Documentation du Centre National de la Recherche 
Scientifique, Paris.) 25s. (1,000 fr.) bound, 22s. 6d. (900 fr.) paper covers. 


There are nine chapters in this book on laboratory methods. The opening chapter 
4 deals with principles of construction and is mainly on the kinematics of apparatus, though 
it also treats the mechanical properties of the commoner materials of construction. The 
* general editor contributes a valuable chapter on vacuum technique, and others are on glass 
| blowing, production of high temperatures, optical technique, light sources and filters, 
| photoelectric cells, recording mechanisms and electrical apparatus. 

t A great deal of information is given in compressed form, and the book will be of interest 
{ to many laboratory workers. From the point of view of readers in this country, it has the 


| incidental interest of teaching us many terms in French which we should find difficult 
(to learn in any other way. ose 


| Ishthara Tests for Colour Blindness, by SuiNopu IsHtHaRa. Reprinted (1948) by 
. meee leewas*&eCo,;: Ltd: ~ 50s; 


Among the different types of colour vision test, the confusion charts are generally 
regarded as the most convenient as a first test. "These charts consist of a series of cards 
| on which a symbol such as a letter or number is formed from a pattern of coloured dots 
} against a differently coloured background. By the correct choice of colours, the numbers 
| may be easily seen by a person with normal colour vision, but may be misread or may be 
/ invisible to a person with defective colour vision, since he may be unable to distinguish the 
colour of the symbol from that of the background. 

A’ successful test not only depends on the correct choice of the confusion colours but 
| also én the shape of the symbols, the size of the dots and so on, since these will affect the 
} ease of recognition of the symbol. The accuracy of reproduction is also vitally important 
| if consistent results are to be obtained from different copies of the same test. One of the 
} most satisfactory of the confusion chart tests, and probably the most widely used, has been 
) the Ishihara test, originating in Japan. During the war there was an acute shortage of this 
test, and an English edition was accordingly published by Messrs. H. K. Lewis in 1943 
and reprinted in 1944. ‘This edition was not quite identical in colouring with the Japanese 
} edition, but in a report on comparative tests carried out in the United States, the English 
} edition was found to be generally satisfactory (7. Opt. Soc. Amer., 1943, 35, 268, 350, 481). 

For some time both the English and Japanese editions have again been unobtainable, 

i but Messrs. Lewis have now issued a third printing (1948) of the English edition. 
| The printing of these charts presents the colour printer with a very difficult problem, but 
comparison between the present printing and that of 1943 shows a very close similarity. 
t Some differences are, however, noticeable, such as a tendency for the colours to be slightly 
stronger in the new printing, although tests on a defective observer did not suggest that 
these differences would have any significant effect on the practical value of the charts. 

An adequate assessment of the charts would require a lengthy investigation with both 
_ normal and defective observers, but the new printing is on the whole so very similar to its 

predecessors that there is little reason to doubt that it should be given the same high rating 
as before. We should like to congratulate the publishers on their enterprise in bringing 
out this reprint, and the printers on their skill in carrying out their difficult task. 

The Physical Society has a special interest in this question, since the Report on 
Defective Colour Vision in Industry prepared by the Colour Group and published by the 
Society recommended the use of the Ishihara test as a general purpose test. At the time 
of publication of the Report, in 1946, it was somewhat unfortunate that the test was 
virtually unobtainable. It is to be hoped that the more extensive testing recommended in 

’ the Report both in schools and in industry may be proceeded with apace now that this new 
copy of the Ishihara test is available. 

Since the term ‘“‘ colour-blindness”’ is now widely regarded as misleading, we should 
like to suggest that any future editions be given the title ‘“‘Ishihara Tests for Colour 
Vision’ or “ Ishihara Tests for Defective Colour Vision ”’, although such a change would 


presumably require the approval of the original publishers, H. M. CARTWRIGHT. 
: W. D. WRIGHT. 
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Frequency Modulation-Volume I, edited by A. N. GoLpsmitH, A. F. Van Dyck 
R. S. Burnap, E. T. Dickry and G. M. K. Baker. Pp. x+515. (RC. 
Review, Radio Corporation of America, R.C.A. Laboratories Divisio 
Princeton, 1948.) $2.50. 


The material of this volume comprises a symposium of 24 articles by various membe 
of the R.C.A. organization, all of which, with one exception, have appeared during the 
period 1936-1947 in either the R.C.A. Review, the Proceedings of the Institute of Ra dio 
Engineers, or various organs of the U.S. radio technical press. "They have been selected — 
so as to cover broadly the main aspects of frequency modulation theory and technique, 
and are arranged according to their subject matter rather than in their chronological order — 
of appearance. After a preliminary group of four papers on the general noise-reducing — 
properties and frequency bandwidth requirements of frequency modulation, the articles 
become more detailed in nature and deal with such subjects as propagation, antennae, — 
‘generation and detection of frequency-modulated waves, and their application to particular 
services such as facsimile. 

In addition to the articles reproduced in full, each section is followed by a number of 
further articles in summary form, and at the end is a bibliography containing some 88 titles 
all from authors connected with the R.C.A. organization. 

The value of this book, which in spite of its title is not to be regarded as a comprehensive © 
text, lies in the assembly into one volume of a number of specialist papers which were 
previously distributed throughout the literature. D. G. REID. 


1 


SMITH, A. F. VaN Dyck, R. S. BurNnap, E. T. DickEy and G. M. K. BAKER 
Pp. x+485. (R.C.A. Review, Radio Corporation of America, Princemmaay 
1948.) $2.50. 


Over the past two decades the Radio Corporation of America has published in al 
R.C.A. Technical Book Series selections of the more important contributions of their staff _ 
to the R.C.A. Review, to the Proceedings of the Institute of Radio Engineers, etc. ‘The 
present book is the eighth volume of this series and the second on the general subject of 
radio, Radio at Ultra-High Frequencies, Volume I, having covered the period 1930-1939. 

The papers in the present book cover rather more than is indicated by the title, and 
include much material which could be classified as v.H.F. and microwave frequency 
technique. ‘The papers are given in seven different sections : antennas and transmission 
lines; propagation; reception; radio relays; microwaves; measurements and components; 
and navigational aids. All the more important papers are published in full, a few are given — 
in abstract, and there are appendices giving a summary of all the papers appearing in~ 
Volume I, and a bibliography covering all the papers published by R.C.A. on U.H.F. from — 
1925 to 1947. ; 


> 
Radio at Ultra-High Frequencies, Vol. II (1940-1947.) Edited by A. N. —_ 


‘ 


The papers given in full have been well chosen and include several important papers on — 
aerial theory and design by G. H. Brown and N. E. Lindenblad and their collaborators. ; 
There is also a very good collection of papers on receivers, including D. O. North’s paper 
on “The Absolute Sensitivity of Radio receivers’’, V. D. Landon’s paper on ‘“‘ The 
Distribution of Amplitude with Time in Fluctuation Noise ”’, as well as several valuable 
papers by E. W. Herold on “ Signal-to-Noise Ratio of U.H.F. Receivers’’. Altogether the 
book may be said to contain a very representative set of papers indicative of many of the | 
latest techniques in the field of U.H.F. radio. 

Many of the recent advances in radio and electronics depend on the research and — 
development work in the higher radio frequencies, and many scientists and engineers who 
have difficulty in gaining access to the original papers will find the collection of papers 
given in this book of real value as the primary source of the contribution of the Radio 
Corporation of America in the field of U.H.F. developments. DENIS TAYLOR. 
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ABSTRACTS FOR SECTION A 


Sie Variation of Resistivity for Mercy and Tin, by E. R. Anpeew. 


RACT. Mezsurcments have been made of the low temperature resistivity of mercury 
z s down 10 6 » diameter and solled tin foils down to 3 p thickness. A contimuous 
ease of resistivity with decrease of diamctes and thickness was found, and it % concluded 
this was duc t9 2 shortening of the mean free path of the conduction electrons by 


ime is cxamaned 5 2 rigorous theory for 2 foil already exists (Fuchs 1933). By comparing 
Z Pagetumcatd cetults with the theory, values of the product Ip are found for both metals; 
1% the mean free path and p the resistivity for the metal in bulk). These values of Ip are in 
facts ee 
n eflect for the two metals. 


Critical Field Measurements on Superrmediactisgy Tox. Fo Foils, by E. R. Anpgew. 
STRACT. The critical fields of superconducting tin foils have been determined from 
measurements in 2 longitudmal magnetic field. Poameranchuk’s relation, 


values of the iength ¥ were found 2s 2 function of temperature ; 20 & the value of ¥ 
estamated to be 12% 10-* cm, rising rapidly close to the transition temperature. Itis 
: aetna endl delanemage a sampoungamy aogier cigeiar ti 
g the validity of Ginsburg’s miterpretation of ¥ as the sum of 4 2nd the surface 
2 sg ae ae wren thee Saeed f ane appuenientcly equal. 


nt the Radiation from Transient Light Sources, by ¥. Boots. 
SSTRACT. A thesetical txeatmnent % given of the problem of analysing the radiation 


E. F. Caldm In the present paper 2 critical examination of this method is made, 
nd exact formuilze desived for special models. 
ei ete ns Oe 2 chee of Secamicoceat partictes, ali at the 


c particle surface, and (b) the particles reflect light diffusely ; (ii) a uniform distribution 
of particles m a2 thin spherical shell 


tasiec collision with the boundary surfaces of the metal. The theory of this effect fora 
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It is shown that the condition for the radiation from a cloud of black-body particles 
to be equal in intensity to that from a black body of the same size and shape is that 
R®ap) <M, where R is the radius of the cloud, a the particle radius, py the particle density 
and M the total mass of the cloud. The cooling of clouds by radiation transfer is also 
examined theoretically. 

The bearing of the results on Caldin’s method is discussed, and more exact methoum 
of analysing the cooling of this form of light source sketched. 


On the Schumann-Runge O, Bands emitted at Atmospheric Pressure, by M. W. 
Feast. 


ABSTRACT. The Schumann-Runge bands of Os, previously analysed in the region 
3100-4200 a. and photographed below 3100 a. under small dispersion, with only a tentative 
analysis, have been obtained from 2000 a. to 5000 a. on sufficient dispersion for the vibra- 
tional analysis to be extended and for the rotational structure to be clearly resolved. 
Conditions for their emission are discussed and the probable process, involving dissocation 
followed by recombination, is outlined. On this basis the results obtained in discharges, 
ozonizers and high-voltage arcs are explained. Flory’s suggestion of predissocation in the 
upper state is reviewed and the experimental evidence found to be against it. The conclu- 
sions reached are applied to considerations regarding the spectrum of the night sky. 


The CN Tail Bands emitted by the Carbon Arc in Air, by M. W. FEast. 


ABSTRACT. Wavelengths and intensities of twenty-one CN tail bands, observed in the 
carbon arc spectrum, are listed and v’, v” values assigned. Nine of the bands have been 
previously reported as occurring in the nitrogen afterglow excitation. The intensities of 
the tail bands are shown to be in accordance with the accepted potential curves for this 
system. 


Ultra-Violet Bands of Naz, by S. P. SINHA. 


ABSTRACT. The ultra-violet bands of Na, between 3600 a. and A 3200 a. have been 
photographed in absorption in the first order of a 21 ft. concave gratiny with a dispersion 
of about 1°3 a/mm. and their vibrational analysis carried out. All the bands belong to a 
single system and can be represented by the equation E 
v= 29342+119-33 (u’)—0-53 (u’)?—159-23 (u”)+0-726 (u”)?+0-0027 (u”)3, 
where u=v+3. . 
The heat of dissociation for the upper state is calculated to be 6510 cm~?. The bands 

are due to transition from the ground !X¢ state to an excited state which is probably 1;¢ 


in character. The dissociation products for the upper state are either 3 °S+4?P or 
32S4+4-32D atoms of sodium. 


On the Equations of Motion of Crystal Dislocations, by F. C. FRANK. 


ABSTRACT. It is shown that when a Burgers screw dislocation moves with velocity v 
it suffers a longitudinal contraction by the factor (1—v?/c?)?, where c is the velocity of 
transverse sound.. ‘The total energy of the moving dislocation is given by the formula 


E:=:Eg/(1 - v/c?)t, where Ey is the potential energy of the dislocation at rest. Taylor — 


disiucations behave in a qualitatively similar manner, complicated by the fact that both ’ 
longitudinal and transverse displacements and sound velocities are involved, 
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